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Preface 



This thesis contains the summary of my research work carried out as a Ph.D. stu- 
dent at the Theoretical Physics Department of the Institute of Nuclear Research of 
the Hungarian Academy of Sciences Debrecen, Hungary and partly at the Theoret- 
ical Physics Institute of the Karl-Franzens Universitat Graz within the framework 
of the fruitful collaboration between the few-body research groups of Debrecen and 
Graz. The new results underlying this thesis have already been presented at in- 
ternational scientific meetings and published in four papers appeared in Journal of 
Mathematical Physics and Physical Review C. [3^, Q. 

As the title "Continued fraction representation of quantum mechanical Green's 
operators" implies my research work is concerned with one of the central concepts 
of quantum mechanical few-body problems. The exploitation of the richness of 
the mathematical theory of continued fractions has enabled us to develop a rather 
general method for evaluating an analytic and readily computable representation 
of Green's operators. This effective representation facilitates the solution of fun- 
damental few-body integral equations. 

Being a theoretical physicist I have always been interested in mathematics. 
Therefore as a second year undergraduate physics student I was enthusiastic to 
accept the research task on two-point Pade approximants offered by my present 
supervisor. Later under his supervision I finished my diploma work on the second 
order Dirac equation. By the time of my PhD studies I got completely infected with 
mathematical physics and my attention was attracted to the topic of continued 
fractions and quantum mechanics. During my work I benefitted a lot from my 
supervisor's knowledge and experience, and I could identify myself with many of his 
thoughts. He showed me how mathematics could be used in practice approaching 
real physical problems. I hope my thesis may serve as an example for how physics 
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has always profited from mathematics. 

Here I would like to take the occasion and thank to everybody who helped me 
during the time I did my research and wrote my thesis. 

Special thanks go to Prof. Zoltan Papp my PhD supervisor for his guidance 
and useful advises and for the excellent atmosphere in which we have been working 
together. 

Some of the work covered by this thesis was done and the results were published 
together with Dr. G. Levai. I thank him for the efforts he invested in our common 
topic. This collaboration was a very instructive experience for me, I could learn 
from him how one could always be optimistic even in the most hopeless moments. 

I am grateful to the Theoretical Physics Department of ATOMKI for providing 
me a peaceful and pleasant working environment. I thank Prof. Borbala Gyarmati, 
who read the manuscript and made numerous helpful suggestions. I am indebted 
to Prof. R. G. Lovas who made the completion of my thesis possible by offering a 
young research fellow scholarship. 

I am also grateful to the few-body group of the Theoretical Physics Institute of 
the Karl-Franzens Universitat Graz, especially to Prof. W. Plessas, for the vivid 
scientific atmosphere. 

Last but not least I wish to thank all the Professors of the Tuesday and Thursday 
4PM open-air seminars for their brilliant lectures on tactics, sport diplomacy and 
football. 
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Chapter 1 



Introduction 



The theoretical description of the microscopic world can be approached following 
two "orthogonal" paths. According to the many-body or field theoretical approach 
the microscopic world is considered as an assembly of many or infinitely many 
interacting objects, where field theoretical or statistical methods can be applied 
in order to describe the system |Q, ^, ||. On the contrary, few-body physicist 
tackle physical systems possessing only few degrees of freedom being consisted of 
a couple of interacting particles and intend to provide a physically complete and 
mathematically well formulated description. 

Few-body systems have played a crucial role in the development of our under- 
standing of microscopic world: atomic, nuclear and particle physics heavily rely 
on few-body models. Nowadays however, few-body problem has an ambiguous 
reputation of being a jungle where non-experts are quickly discouraged and spe- 
cialist enjoy endless debates on technical improvements concerning equations and 
mathematical physics issues. Furthermore much attention has been paid only to 
computational consequences and little concern about the underlying physics. This 
is of course a false impression which originates from the non-trivial nature of the 
problems studied. The goal of the few-body physics community, namely giving a 
complete and mathematically correct description of few-body systems, automaticly 
requires the consideration of mathematical issues, hence this field offers an excel- 
lent playground for mathematical physics. This thesis, following the above ideas, 
hopes to contribute to the magnificent results achieved by the few-body approach. 
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CHAPTER 1. INTRODUCTION 



Microscopic few-body physics certainly was born together with quantum me- 
chanics in the pioneering work of Bohr, Heisenberg and Schrodinger trying to de- 
scribe simple quantum systems, like the Hydrogen atom. The field initially de- 
veloped as a part of nuclear physics as the title of the first few-body conference 
(Nuclear forces and the few-nucleon problem |Q held in London 1959) suggests. 
Since then few-body physics has expanded to incorporate atomic, molecular and 
quark systems. For example, the two-electron atom had already been attacked 
by Hylleraas |j| in 1929 and conquered by Pekeris Q only in 1959. The theoret- 
ical foundations of few-body quantum physics were laid down by Lippmann and 
Schwinger Gell-mann and Goldberger || by developing formal scattering theory 
of two-particle systems. The first attempts to extend the results to multichannel 
processes and more than two particles led to unsound mathematical formalism and 
non-unique solutions. The rigorous theory of few-body systems was given by Fad- 
deev |^| who proposed a set of coupled integral equations which have a unique so- 
lution for the three-body problem. Having the correct few-body theory much effort 
has been invested into the development of numerical methods. The first numerical 
solution of the Faddeev equations for three spinless particles interacting with local 
potentials was achieved by Humberston at al. [|l(| in 1968. Since then a lot has 
been achieved due to the unbelievable development in computational power and 
the several extremely effective new methods having been developed. Among others, 
configuration and momentum space Faddeev calculations jn], [l^, |l]|, the hyper- 
spherical harmonics expansion method Jli), the quantum Monte Carlo method 
[HH , the Coulomb-Sturmian discrete space Faddeev approach Jl(| and the stochas- 
tic variational method |l^| have been used to study few-body bound, resonant and 
scattering phenomena with great success. 

The fundamental equations governing the dynamics of few-body physical sys- 
tems, like the Lippmann-Schwinger equation and the Faddeev equations, are for- 
mulated in terms of integral equations. Integral equation formalism has a great 
advantage over the equivalent traditional differential equations because of the few- 
body boundary conditions are automatically incorporated into the integral equa- 
tions. This is the reason why integral equation methods are in favour when scatter- 
ing problems with complicated asymptotic behaviour are considered. Nevertheless 
for cases, where the asymptotics of the wave function is well known, differential 
equation approach can perform outstandingly (see i.e. the stochastic variational 
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method [§])). 

In spite of the fundamental merits of integral equations, their use in prac- 
tical calculations is usually avoided and various approximations to the few-body 
Schrodinger equation are preferred instead. The reason for this is certainly that the 
integral equations contain not the usual Hamiltonian operator, but its resolvent, the 
Green's operator in their kernel. The evaluation of the Green's operator is much 
more complicated than the direct treatment of the Hamiltonian using standard 
tools of theoretical and mathematical physics. The determination of the Green's 
operator is equivalent to the solution of the problem characterized by the corre- 
sponding Hamiltonian. In standard textbooks on quantum mechanics jl^] Green's 
operators are introduced on the level of fundamental equations and their general 
properties are extensively studied. Formal scattering theory |l9|, [20|, thus few-body 
quantum mechanics can be formulated upon the notion of Green's operator. The 
Green's operator, similarly to the wave function, carries all the information of the 
physical system. Consequently the determination of the Green's operator, as the 
central concept of few-body quantum mechanics, is of extreme importance. 

The suitable choice of the Hilbert space representation can facilitate the de- 
termination of the Green's operator. In this respect the momentum space rep- 
resentation is rather appealing as the free Green's operator is extremely simple 
there. This is the main reason why momentum space techniques are so frequently 
used and also why they are capable of coping with complicated integral equations 
(for a review see [^||). On the other hand the use of discrete Hilbert space basis 
representation is often very advantageous because it transforms the integral equa- 
tions into matrix equations. The harmonic oscillator (HO) functions and the 
Coulomb-Sturmian (CS) functions are good examples of discrete Hilbert space ba- 
sises. The free Green's operator can also be given analytically between harmonic 
oscillator states |^2|. This allowed the construction of a flexible method |23] , p4| 
for solving the Lippmann-Schwinger equation, which contains the free Green's op- 
erator in its kernel and provides a solution with correct free asymptotics in HO 
space. The CS basis representation of the Coulomb Green's operator in terms of 
well computable special functions was derived by Papp in Ref. |^5| , where he could 
perform complicated analytic integrals of the Coulomb Green's function and the 
Coulomb-Sturmian functions making use of the results of Ref. [^| . 

The CS representation of the Coulomb Green's operator forms the basis of 
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a quantum mechanical approximation method developed recently J2j], |2?j, for 
describing Coulombic systems. The strength of the method is that the Coulomb- 
like interactions in the two-body calculations are treated asymptotically correctly, 
since the Coulomb Green's operator is calculated analytically and only the asymp- 
totically irrelevant short-range interaction is approximated. This way the correct 
Coulomb asymptotics is guaranteed. The corresponding computer codes for solving 
two-body bound, resonant and scattering state problems were also published p9[ . 
The method has been extended to solve the three-body Coulomb problem in the 
Faddeev approach. In this formulation of the Faddeev equations the most crucial 
point is to calculate the resolvent of the sum of two independent, thus commuting 
two-body Coulombic Hamiltonian. This resolvent is given as a convolution integral 
of two-body Green's operators. Therefore the evaluation of the contour integral 
requires the analytic knowledge of the two-body Green's operators. So far good 
results have been obtained by Papp for bound state and for below-breakup 
scattering state three-body problems [ |30| , 

In practice there is no general procedure for how to determine Hilbert space 
representation of a Green's operator. All the Green's operators mentioned before 
require separate and detailed investigation sometimes based on very specific con- 
siderations. In a recent publication |3lJ we have proposed a rather general and 
easy-to-apply method for calculating discrete Hilbert space basis representation of 
Green's operators belonging to some class of Hamiltonians. We have shown that 
if in some basis representation the Hamiltonian possesses an infinite symmetric 
tridiagonal (Jacobi) matrix structure, the corresponding Green's operator can be 
given in terms of a continued fraction. The procedure necessitates the analytic 
calculation of the matrix elements of the Hamiltonian as the only input in the 
analytic evaluation of the Green's matrix in terms of a convergent continued frac- 
tion. This method simplifies the determination of Green's matrices considerably 
and the representation via continued fraction provides a readily computable and 
numerically stable Green's matrix. The Green's operator of the non-relativistic 
two-body Coulomb problem and the D-dimensional harmonic oscillator problem 
was evaluated utilizing this method in Ref. J3l]]. An exactly solvable potential 
problem, which provides a smooth transition between the Coulomb and the har- 
monic oscillator problem, was considered in Ref. |Q. A suitable basis was defined 
in which the Hamiltonian of the potential problem appeared in tridiagonal form. 
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With the help of the analytically calculated matrix elements of the Hamiltonian 
the determination of the Green's matrix was straightforward. The Hamiltonian of 
the radial Coulomb Klein-Gordon and second order Dirac equation was shown to 
possess infinite symmetric tridiagonal structure in the relativistic CS basis. This 
again allowed us to give an analytic representation of the corresponding relativistic 
Coulomb Green's operators in terms of continued fraction [Q. The continued frac- 
tion representation of the Coulomb Green's operator was utilized to give a unified 
description of bound, resonant and scattering states of a model nuclear potential 
using a quantum mechanical approximation method |^3| devised to solve the 
Lippmann-Schwinger equation. 

The layout of this thesis is the following: 
This introductory chapter is followed by a chapter devoted to the concept of the 
Green's operator. Here few-body quantum theory is built up around the Green's 
operator. In Chapter || a review of tridiagonal matrices, three-term recurrence rela- 
tions and continued fractions is given. The method of calculating continued fraction 
representation of certain class of Green's operators is presented in Chapter ||. This 
chapter also covers the calculation of continued fraction representation of the D- 
dimensional Coulomb, D-dimensional harmonic oscillator, the generalized Coulomb 
and the relativistic Coulomb Green's operators. In the last chapter two typical ap- 
plications of the Green's operator in few-body physics are delivered. First the 
continued fraction representation of Coulomb Green's operator is used to present 
a unified description of bound, resonant and scattering states of a model nuclear 
potential demonstrating the central role of Green's operators in few-body quantum 
mechanics. Then the binding energy of the Helium atom is calculated by solving 
the Faddeev-Merkuriev equations for the atomic three-body problem in Coulomb- 
Sturmian Hilbert space representation. Chapter |5] and [| contain the new results of 
the author. Finally, summary and bibliography closes the thesis. 



Chapter 2 



The quantum mechanical 
Green's operator 



Expressions like Green's operator, Green's function and propagator are frequently 
used in different fields of physics and mathematics sometimes with confusing mean- 
ings. For example the many-body Green's function G a p(xt,x'tf) studied in quan- 
tum many-body problem is defined as the ground state expectation value of fields 
operators in the second quantized formalism Whereas the propagator of rel- 
ativistic quantum field theory is often referred to as the Green's function as well 



The Green's operator of few-body quantum physics corresponds to the definition 



2.1 Definition 

Suppose H is a Hermitian differential operator in some Hilbert space. The G{z) 
Green's operator is defined as the resolvent operator of H 



I 



of mathematical physics 
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G(z)(z -H) = (z- H)G{z) = 1 



for z £ C/a(H), 



(2.1) 
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2.2. BASIC PROPERTIES OF GREEN'S OPERATOR 



where a(H) denotes the spectrum of the Hermitian operator H and 1 is the identity 
operator. So, the G(z) = (z — H)^ 1 Green's operator is defined on the whole 
complex plane C except for the spectrum of H . 

The {\r)} coordinate basis representation of the (O) operator equation 



S(r - r') = (r\I\r'), J dr|r)(r'| = 1 (2.2) 

G(r, r'; z) = (r\G(z)\r'), 5(r - r')H(r) = (r\H\r') (2.3) 

leads to an inhomogeneous differential equation 

[z-H(r)]G(r,r';z) = 6(r-r') (2.4) 

for G(r,r';z). In mathematical physics the Green's function corresponding to the 
linear Hermitian differential operator H(r), where r S f2, and to the complex 



variable z is defined as the solution of the inhomogeneous equation (2,4) subject 
to certain boundary conditions for r, r 1 on the surface of the domain f2. In fact the 
Green's function is required to satisfy the same boundary condition as the wave 
function. Consequently the G(r, r'; z) Green's function of mathematical physics can 
be viewed as the coordinate space representation of the G(z) resolvent operator. 

In what follows, under the notion of few-body Green's operator the resol- 
vent of the Hamiltonian is to be understood. In the next chapters Hilbert space 
basis representation of Green's operators will be investigated. 



2.2 Basic properties of Green's operator 

The determination of the Green's operator represents the solution of the underlying 
differential equation. Therefore the evaluation of the Green's operator correspond- 
ing to the Hamiltonian H is equivalent to the complete description of the system 
characterized by H. The Green's operator carries all the information about the 
physical system. Once some representation of G(z) is known one can gain the 
eigenvalues and the eigenfunctions of the Hamiltonian, construct projection opera- 
tors or determine the density of states, etc,. The time evaluation operator can also 
be obtained by performing a Fourier transform of the Green's operator. 

Let {|</>n)} and {a n } denote the complete set of orthonormal eigenstate and the 
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corresponding eigenvalues of the Hermitian operator H 

ff|0n) = a n \4> n ) . (2.5) 

Utilizing the completeness of the eigenstates, the Green's operator can be written 

as 

G[z) = ^MM + [JMM, (2 . 6) 

n 

where J2 n anc ^ / ^ n indicates summation over the discrete and the continuous 
spectrum, respectively. 

Since H is Hermitian, all of its a n eigenvalues are real, which means that G(z) 
is an analytic function in the complex z-plane except at those points and intervals 
of the real axis which corresponds to the spectrum of H. According to ( ^lj| ) G{z) 
possesses simple poles at the positions of the discrete eigenvalues of H, hence the 
poles of G(z) give the discrete energy eigenvalues of the Hamiltonian. The residue 
at the a n pole provides information on the corresponding non-degenerate \(j> n ) 
eigenvector according to the projection operator 

\4>n){4>n\ = ^ i j dzG{z), (2.7) 



here the contour encircles the a n and only the a n eigenvalue. Equation (2.7) can 
be used to construct the wave function. 

If A e K belongs to the continuous part of the spectrum, then G(A) is usually 
defined by a non-unique limiting procedure 

G+(A) := lim G(X + ie) (2.8) 

£->0+ 

G~(X) := lim G(X-ie) (2.9) 

because of the existence of the limits. Thus the continuum spectrum of H ap- 
pears as a singular line (branch cut) of G(z). Due to the self-adjointness of the 
Hamiltonian the Green's operator exhibits the important property 



G{z*) = G(z)* 



(2.10) 



2.3. THE GREEN'S OPERATOR IN SCATTERING THEORY 
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Therefore G + (A) and G (A) are connected by 

G-(A) = (G+(A))% (2.11) 

from which we have 



ReG-(A) = ReG+(A), ImG"(A) = -ImG+(A) . (2.12) 

The knowledge of ReG ± (A) allows us to obtain the density of states at A 

A(A) = T-HrrfG^)]} . (2.13) 

The time development of the \i/>(to)) state belonging to the Hamiltonian H at 
t = to can be obtained as a solution of the time dependent Schrodinger equation, 
and it is governed by the U(t,to) time evolution operator according to 

\m) = u(t,t )m )), (2.i4) 



where the time evolution operator 

U(t,t ) = e -*^(*-*o) (2.15) 



can be expressed in terms of the Green's operator as 
[/(Mo) = ih 



OO 1 



2ir 



(2.16) 



Equation ( 2.1 6| ) shows that the U (t, to) propagator can be determined as the Fourier 
transform of the G ± Green's operators. 



2.3 The Green's operator in scattering theory 

The formulation of quantum scattering theory is heavily relied on the concept 
of Green's operator. It is concerned with the mathematical description of the 
scattering process during which the incident beam of particles, usually prepared by 
some accelerator, enters into interaction in the scattering region then the scattered 



10 



CHAPTER 2. THE QUANTUM MECHANICAL GREEN'S OPERATOR 



beam is observed by detectors located at large distance from the interaction area. 

Far from the scattering region the incoming and the scattered beam can be 
thought of as bunch of particles propagating freely, therefore they can be charac- 
terized by the \ipi n ), IVWt) asymptotic states whose time evolution are governed 
by the H° asymptotic Hamiltonian as 

\An/out(t)) =Uo(t)\An/out(to)), U (t) = e" '* * '* . (2.17) 

On the other hand the time evolution of the scattering process is fully determined 
by the H = H + V total Hamiltonian involving the scattering potential as well: 

m))=U(t)m Q )), U(t) = e-i m , (2.18) 

here \i/>(t)} denotes the actual scattering state vector. However, during scatter- 
ing experiments one can only measure the \ipi n ) incoming and the \tpout) outgoing 
asymptotes of \ip). Therefore in scattering theory all information about the scat- 
tering process should be extracted from the asymptotic behavior of \ip). 

At t = let (\x— )) denotes the \i/)(t)) scattering state which evolved 

from the \ipi n ) = \4>) incoming asymptote (evolved to the \ipout) — \x) outgoing 
asymptote). The \<f>+) and \x~ ) scattering states are related to their asymptotes 
by the equations 

= n+M (2.i9) 

lx-) = «-lx>, 

where the fl± M0ller operators are defined as the limits of the time evolution 
operators 

Q± = lim U(tfll a (t). (2.20) 

t — >=foo 

The existence of the M0ller operator can be proved if the scattering potential 
satisfies certain conditions (i.e. not too singular at the origin, sufficiently short 
ranged and smooth, see [Q). The main goal of the scattering description is to relate 
the outgoing asymptote to the incoming one without reference to the experimentally 
unknown actual orbit. Utilizing the f2LJ7_ = 1 isometric property of the M0ller 
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operator from Eq. ( 2.19 ) we obtain 

|Vw) = nln+lVfo.) := S\ip in ), (2.21) 

here S := is the scattering operator, which carries all the experimentally 

important information about the scattering. Once the operator S is determined 
the scattering problem is solved. 

In order to determine the S operator it is convenient to use the { \p) } momentum 
representation. It is important to note here that the \p) plane wave momentum 
eigenstates do not represent physically realizable states, in fact they are only used 
in the expansion of physical states of the scattering process as 



= / d'p V(p)|p). (2.22) 
Similarly the scattering state |0+) can be expressed as 

= Q + \d>) = [ d 3 p <f>(m+\P) = I d 3 P <t>(p)\P+), (2-23) 



where the \p+) = improper (i.e. not normalizable) vector is called the 

stationary scattering vector corresponding to the \p) incoming plane wave. 

In momentum representation a careful mathematical analysis leads to the very 
important result 



a+\<f,) = |0+) = lim tf(t)t[T o (t)|0} 

t — > — OO 

p — OO 

= \4>)+i dr U{T)^VU a {T)\4>) 



-oo 



(2.24) 

lim / dr e +eT U(ryVU Q {T)\<t>) 



= \4)+i lim / d> G{E P + ie)V\p) (p\4>) 

e-> 0+ J 

which establishes the connection between the Q+ M0ller and the G + {z) Green's 



operator. After expanding \4>+) and \<f>) in Eq. (2.24) in terms of the improper 
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states \p) and \p+) we obtain for the stationary scattering vector \p±) 



Making use of the 



\p±) = |p) + lim G(E p ± ie)V\p) 

e— >0+ 



A" 1 = B- 1 + B~ r (B - A)A~ 1 



(2.25) 



(2.26) 



operator identity with A = z — H and B = z — H°, then with their interchange, 
we can derive the 



G(z) = G°(z) + G°{z)VG(z) 
G(z) = G°(z) + G(z)VG°(z) 



(2.27) 



resolvent equations, which relate the G° and the G operators. With these equations 



Eq. (2.25) can be recast in the form 



\p±) = \p)+ lim G°(E p ± ie)V\p±). 



(2.28) 



This type of integral equations were first formulated by Lippmann and Schwinger 
[Q, consequently Eq. ( 2.27 ) and Eq. ( 2.28 ) are known as the Lippmann-Schwinger 
equations for G and \p±) , respectively. According to ( 2.25|) and ( 2.28| ) the evalua- 
tion of the stationary scattering vectors necessitates the determination of the full 
Green's operator or the solution of the ( [2.28 ) Lippmann-Schwinger equation. 



Being able to calculate |p±), one can obtain any scattering information. For 
example, the (p'\S\p) scattering matrix can be written as 

(p'\S\p) =S 3 (p'-p)- 2mS(E p> - E p ) t(f ^ p). (2.29) 

Here t(p ' <— p) can be calculated using the |p±) scattering states 

tip' = (p'\V\p+) = (p'- \V\p), (2.30) 
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or using the Green's operator 

tip' «- V) = (p'\V + VGV\p). (2.31) 

The (V + VGV) operator is called the T operator. Its tip ' <— p) matrix element is 
referred to as the on-shell T matrix, because in Eq. ( |2.2g| ) it appears together with 
the S(E p i — E p ) delta function, and thus defined only for the p' 2 = p 2 "on-shell" 
energy case. 

The experimentally off important differential scattering cross section da/dfl can 



be derived from the ( |2.26| ) formula for the scattering matrix, and has the form 

^-l/V-JOl*. (2.32) 

where f(p' <— p) is the scattering amplitude, which is just a trivial factor times 
the on-shell T matrix 

f{p' «- p) = -(2n) 2 m t(p' <- p). (2.33) 



Finally we quote here the asymptotic form of the coordinate representation of 
the stationary scattering vector 



lim (r\p+) = (27r)" 



e lpr + f(pf^p) 

r 



(2.34) 



which reflects the traditional definition of scattering states as the sum of an incident 
plane wave and a spherically spreading scattered wave, whose f(pr <— p) amplitude 



determines the scattering cross section via (2.32) 



In this short review it was demonstrated that the knowledge of the Green's 
operator enables one to determine the scattering states, which is equivalent to 
the solution of the scattering problem and makes possible the calculation of any 
scattering quantities. 
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2.4 Dunford-Taylor integrals of Green's operators 

From the theory of linear operators |3(| we know that the analytic function of 
a bounded self-adjoint operator A can be defined as the Dunford-Taylor integral 
representation 

f(A) = ^ I dzf(z)(z-A)-\ (2.35) 

2VT1 JC 

where C encircles the spectrum of A in positive direction and / should be analytic 
on the domain encircled by C. This definition can be considered as the operator 
equivalent of the Cauchy integral formula. 

In few-body physics one often encounters the problem of determining Green's 
operators corresponding to Hamiltonians composed of two independent operators 
with known resolvents. This is the typical situation when the Green's operator of 
separable subsystems are known (e.g. think of a three-dimensional motion decom- 
posed to the sum of a two- and a one-dimensional motion). The question, whether 
it is possible to determine the Green's operator of the composite system making 
use of the sub-Green's operators naturally emerges. The answer can be given by 



utilizing the (2.35) Dunford-Taylor integral representation. 

Let hi,li2 be two independent (hence commuting) bounded linear operators 
with the analytically known gi(z) = (z — hi)^ 1 , g 2 {z) = (z — resolvents. 



Then applying the ( 2.35 ) definition for /(A) = (z — hi — A) -1 and A = h%, the 
G{z) = (z — H)^ 1 Green's operator of the H = hi + h 2 composite Hamiltonian can 
be evaluated by performing a convolution integral of the gi(z), gi{z) sub-Green's 
operators 

G(z) = ^ I dz' 9l (z-z') g 2 (z'), (2.36) 
J c 

where the contour C should encircle, in counterclockwise direction, the spectrum 
(discrete and continuous) of /12 without penetrating into the spectrum of h\. In Fig. 



2.1a typical integration scenario in the complex z'-plane is shown for a bound state 



system (z < 0). We note that among others, Bianchi and Favella |3j proposed a 



convolution integral of the type of (2.3f) for determining the resolvent of the sum 
of two commuting Hamiltonian. 
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The ( 2.35 ) integal representation can also be used to construct other well-known 
operators as a Dunford-Taylor integral of the Green's operator. For example, the 
identity operator 

1 = — * dz (z-h)-\ (2.37) 



the Hamilton operator 



h= —. f dz z(z- h)' 1 , (2.38) 
J c 

or the time evolution operator 

exp(-ifti) = — * dz exp(-izt) (z - hy 1 (2.39) 
2tti J c 

can be given as analytic integrals of the g = (z — h)^ 1 Green's operator. 
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complex z' - plane 




— < 


spectrum of g^z-z') 


/ spectrum 


of g 2 (z') 


Z ' 


XXX 


► 

contour C 



Figure 2.1: The integration path C in the complex z'-plane for the convolution 
integral of the ( |2.36; ) Green's operator. 



Chapter 3 



Tridiagonal matrices, recurrence 
relations and continued fractions 

Based on the mathematical literature |?9[ [ll], this chapter is dedicated to 
the survey of the mathematical machinery being used in the forthcoming part of 
the thesis. First a truncated inverse of an infinite symmetric tridiagonal matrix is 
presented. Tridiagonal matrices inherently imply three-term recurrence relations 
for the matrix elements of their inverse. For this reason the basic properties of three- 
term recurrence relations are outlined. The final section discusses the fundamentals 
of the theory of continued fractions. The analytic continuation and an important 
type of transformation of continued fractions is investigated. The intimate relation 
of continued fractions and three-term recurrence relations is revealed by Pincherle's 
theorem. 

3.1 Infinite tridiagonal matrices 

Infinite rank matrices naturally emerge in quantum mechanics when physical ob- 
servables are represented on a discrete infinite dimensional basis. Moreover the 
matrix representation of many physical operators are tridiagonal in a certain basis. 
In fact, some computational methods, like the Lanczos method |39|, are based 
on the scheme of creating a basis that renders a given Hamiltonian tridiagonal. 
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However we are always limited to consider finite rank truncations of infinite matri- 
ces. Here we derive a formula [[[l) for the n-th leading submatrix of the inverse of 
an infinite symmetric tridiagonal matrix. The underlying idea originally came up 
in the work of Pal |Q| , and we give here the proof. 

Let A := [oy], i,j = 1, 2, . . . , oo, be an infinite symmetric tridiagonal (Jacobi) 
matrix: 

a«_i=/3j_i, a,i=aj, a u+1 — fa, a tj = for |i-j'|>l. (3.1) 
Let B := [6y], i,j = 1,2, . . . , oo denote the infinite rank matrix, for which 

oo 

1 = AB = BA, % = a *kb k] ■ (3.2) 
fe=i 

Then the following theorem can be stated. 



Theorem 1. Let us define the n-th leading submatrices of A and of B as 

An := [dij], B n := i,j = 1,2,... ,n, respectively. Then 



B„ = {A n + F n ) 



(3.3) 



where F n := [fij], i, j — 1, . . . , n, is given by 



fij = S 







pln+1 



i,3 <n 
i = n,j < n 
j = n,i < n 
i=j = n 



(3.4) 



Proof of Theorem Utilizing the ( |3.l| ) tridiagonal property of A, the ( |3.2| ) definin- 
ing relation of B for i, j < n can be written as 



5ij — a ikbkj + SinUnn+lbn+lj , hj < n , 



(3.5) 



k=l 
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which we can rewrite in the form 



n / h \ n 

Sij — ( dikbkj + Sinain+lbn+ljSknJ-^- ) = 

fe=l V ^ / fc=l 



aik + SinSk n ain+i ) bjy. 

(3.6) 

The inverse of a symmetric tridiagonal matrix possesses the following property |Q : 

^={ Pi ^ * \-J. ■ (3-7) 
I PjQi, if 3 < i 



Therefore 

b n +lj _ PjQn+1 _ Plln+1 _ 
bnj Pjqn Piq n bin 

and follows the statement of the theorem: 



(3.8) 



% = ( a,ik + faL«m+i ) b kj = ^2{aik + fik)bkj- (3.9) 



□ 



Finally we give a schematic picture of the theorem. The infinite symmetric 
tridiagonal matrix A and its n-th leading submatrix A n can be illustrated as 



.4 



( A n ■ \ 



V 



A, 



(3.10) 



/ OCXG 

Similarly the infinite B and the truncated B n general matrices can be plotted 
( B n ■ ■ • \ ( ■ ■ ■ -\ 



B = 



B n = 



V / oox 



(3.11) 
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With the above schematic representation of matrices, the ( pT^ ) statement of The- 
orem |l| can be recast in the form 



/. . . .\ 



V • • • • J 



(■ ■ 



\ 



■ ■ ) 



( 



\ 



fn,n J 



(3.12) 



where /„„ = a nn+1 b ln+1 /b ln . 



According to Eq. (3.12) the truncated inverz B n is readily calculable from the 



dij tridiagonal matrix elements, provided the &i n +i/&i n quotient is at our disposal. 



3.2 Three-term recurrence relations 



Due to Eq. ( ]3.2| ) the matrix elements bij of the inverz of an infinite symmetric 
tridiagonal matrix automatically satisfy a three-term recurrence relation in both 
indexes 

5ij = Pi-ih-ij + o^ibij + Pih+ij, for i = 1, . . . , oo, j fixed, (3.13) 
8ij = Pj-xhj-i + ctjbtj + ffjbij+i, for j = 1, . . . , oo, i fixed. (3.14) 

Let us consider the (homogeneous) three-term recurrence relation 

x n +i = b n x n + a n x n -i, a n ^ 0, n= 1,2,3,..., (3.15) 

where a n ,b n complex numbers are the coefficients of the recurrence relation. A 
sequence of complex numbers {i„}™ is called a solution of ( |3.15| ) if its x n elements 



satisfy the equation for all n £ N. The set of all solutions {a;„}^L of ( 3.15 ) form 
a linear vector space of dimension two over the field of complex numbers, and the 
zero vector of this space is the {0}^L trivial solution. A distinguished type of 
solution is the minimal solution. If there exists a non-trivial solution {m„}™ =0 and 



another solution {d n }^L of (3.15) such that 



(3.16) 
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then {m n } is called the minimal (subdominal) solution. A solution which is not 
minimal, like {d n }, is called the dominant solution. In general, a recurrence relation 
may or may not have a minimal solution, and it is clear that if {m n } is a minimal 
solution then {m 1 } :— c{m} for Vc 7^ is a minimal solution as well. Moreover any 
other non-trivial {y n } solutions are dominant since they can be written 

{Vn} = ci{m„} + c 2 {d„}, c 2 ^ 0, (3.17) 

where {d n } is a dominant solution. The existence of a minimal solution is strongly 



related to the convergence of a continued fraction (see Section 3.3) 



Suppose the x and x\ elements of the ( 3.15| ) three-term recurrence relation 



are at our disposal and we want to apply the recurrence relation in a direct way 
to determine {x n }^ =0 recursively. However, for a minimal solution this method 
does fail in practice. Gautschi |If| pointed out that a direct recursive calculation 
of a minimal solution of a three-term recurrence relation is numerically unstable. 
A surprising demonstration of this instability can be found in |4^| on page 219. 

We note here that three-term linear recurrence relations are closely linked to 
linear difference equations of order two. 



3.3 Continued fractions 



The study of continued fractions, i.e. mathematical expressions of the form 



ai 

b Q + (3.18) 

£12 

61 + 
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started in the 16th century with a work of Bombelli |4(| in 1572, who wrote down 
the first finite continued fraction 

4 

3 + =3.6 (3.19) 



to approximate \/l3 w 3.6055. The traditional form ( 3.1§| ) of a continued fraction 
is often written more economically as 

, ai a 2 a 3 a n 
h(z) + — — — —. 3 " 2 °) 

Oi + b 2 + b 3 + +b n 

The first continued fractions were constructed from integers and were used to ap- 
proximate various algebraic numbers, among others the n. From the 17th century 
on the continued fraction evolved into an important tool of number theory. First 
Schwenter, Huygens, Wallis latter Lagrange, Legendre, Gauss and Galois made 
extensive use of regular continued fractions, continued fractions with integer coef- 
ficients, in their number theoretical studies. 

Euler extended the notion of a continued fraction by generalizing its coefficients 
as functions of complex variables rather than simple numbers. This way continued 
fraction expansions could be used as special tools in the analytic approximations 
of special classes of analytic functions. Following Euler's results a long list of bril- 
liant mathematicians from the 19th century, names like Laplace, Jacobi, Riemann, 
Stieltjes, Frobenius and Tchebycheff made flourish the analytic theory of continued 
fractions within the framework of the classical 19th century analysis. 

Mathematical physics had also discovered the continued fractions. Continued 
fraction expansion of special functions, solution of three-term recurrence relations 
and Pade approximants represent important field of applications. The available 
continuously increasing computational power has also contributed to the present- 
day success of continued fractions in mathematical physics. 

Rigorously a continued fraction is defined as an ordered pair of the type 

(({M*)}, {&»(*)}), {/»(*)}), (3-21) 
where {a n (z)}f and {frn(z)}o°, with all a n (z) ^ 0, are two sequences of complex 



3.3. CONTINUED FRACTIONS 



23 



valued functions defined on the region D of the complex plane. The complex values 
of a n (z) and b n (z) are called the n-th partial numerator and partial denominator of 
the continued fraction, or simply the coefficients or elements. The {fn(z)} sequence 
of complex functions is given by 

fn(z) = S n (0,z), n = 0,1,2,..., (3.22) 



where 

S n (w n ,z) = S n -l(s n (w n ,z),z), So (wo , z) = Sq (wq , z) , (3.23) 



with the 

Sn(w n , z) — ""(^ — n > 1, s (w 0l z) — b (z) + wq, (3.24) 
b n {z) + w n 

linear fractional transformation. 

Here S n (w n ,z) is called the n-th approximant of the continued fraction with 
respect to the {w n }^L complex series. Applying recursively the ( |3.24 ) linear 
fractional transformation the n-th approximant of the continued fraction can be 
written as 

Sn(w n ,z) = b (z) + — -- —— ■■■ . 3.25) 

b 1 {z) + b 2 {z)+ +b n (z) + w n 

Similarly, for the S n (0, z) approximant we got 

S n (0, z) = b Q (z) + K ( $>) = b (z) + + ■ ■ ■ |f| , (3.26) 

i=i\bi(z)J b 1 {z) + b2(z)+ +b n (z) 

where the new notation K^ =1 (ai/bi) was introduced for the approximant. Subse- 
quently for the (({a n (z)}, {b n (z)}), {f n {z)}) continued fraction one of the 

, . , a\(z) clq,{z) , , . oo / a n (z)\ 

W + rH,TTT , ••• . or & o^ + ^ hrH 3 - 27 

6i(z) + 6 2 (z) + n=l\ v 0n(2 : )/ 

notations will be used. 



The convergence of a continued fraction means the convergence of the sequence 
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of approximants S n (w n , z) to an extended complex number 

f(z) = lim ) . (3.28) 

n — >oo 

If f(z) exists for two different sequences of {w n } then f(z) is unique. For a detailed 
discussion of convergence results see Chapter IV of |Q . 

There are several algorithm to compute the n-th approximant of a continued 
fraction. The backward recurrence algorithm consists of summing-up the 



fraction starting at its tail. This method is proved to be numerically stable. The 
main drawback of the backward recurrence algorithm is that it needs to be re- 
calculated for each approximant. An alternative method, the forward recurrence 
algorithm is based upon the following 

B n (z) + B n -i{z)w 

representation of S n (w n , z), where A n (z), B n (z) are called the n-th numerator, de- 
nominator, respectively. The n-th numerator and denominator satisfy a recurrence 
relation namely 

AJz) = b n (z) + A n _i(z) + a n (z)A n - 2 (z) (3.30) 
B n (z) = b n (z) + B n _-L(z) + a n (z)B n - 2 (z), (3.31) 

with A-i = 1,Aq = &o,-B-i = 0, Bq = 1. In contrast with the forward algorithm, 
one can easily obtain S n +i(w n +i, z) from S n (w n , z) using the backward algorithm, 
but on the other hand, this method is less stable than the first one. 

A special class of continued fractions for which the limits 

lim a n {z) = a(z) and lim b n (z) = b(z) (3.32) 

n — >oo n — » oo 

exist for all Z € D is called the limit 1-periodic continued fraction. The more 
general case is the limit-k periodic continued fraction for which we have 

lim akn+ P {z) = a p (z) and lim bkn+ P (z) = b p (z) for p = 1, 2, . . . , k . 

n — >oo n — >oc 

(3.33) 
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Limit periodic continued fractions play an important role in the analytic theory of 
continued fractions, since most of the continued fraction representation of special 
functions are limit periodic continued fractions. The convergence properties of 
limit periodic continued fractions are determined by the behavior of their tail by 
means of the w±(z) fixed points of the limit linear fractional transformation 

s(w,z) = lim s n (w n ,z) — ^ - — . (3.34) 

n— >oo o(z) + W 

The w±(z) fixed points are given as the solutions of the s(w) — w fixed point 
equation 

w= tt4t => «"±0O = "K«)/2 ± >/(6(«)/2) a + a(z). (3.35) 
o(z ) + w 

The fixed point with smaller modulus is called attractive fixed point, while the 
other one is named as the repulsive one. Since w±(z) represent the tail of a 
limit 1-periodic continued fraction we can accelerate the convergence using the 
attractive fixed point in the approximant S n (w,z) Et|]. 



3.3.1 Analytic continuation of continued fractions 

The idea of analytic continuation of a continued fraction bo(z) + K(a n (z)/b n (z)) 
by means of an appropriate choice of w n (z) for the S n (w n {z)) approximant was 
proposed by Waadeland [Q and later recalled by Masson . By examining limit 
periodic continued fractions they sought a modification of the tail of a continued 
fraction which led to its analytic extension. 

If a continued fraction converges in a certain complex region D, then in many 
cases it is possible to extend the region of convergence to a larger domain D* D D, 
where D* depends on the choice of the functions w n {z). 

In the case of limit 1-periodic continued fractions the analytic continuation 



Jd* (z) of the continued fraction f(z) in ( |3.2g| ) is defined with the help of the fixed 



points w±(z) of Eq. ( |3.35| ) as 

Jd'{z) = lim S n {w±(z),z). (3.36) 
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In Eq. ( 3.36 ) the analytic expressions for the w±(z) fixed points are continued 
analytically and then employed to sum up the continued fraction, this way providing 
the analytic continuation of the continued fraction. 



Bauer Muir transformation 

The numerical computation of approximants S n (w±(z), z) might be unstable, spe- 
cially for z belonging to the extended region, which leads to unsatisfactory con- 
vergence. This problem can be overcome by using the Bauer-Muir transformation 
of a continued fraction (see eq. |4lj) We note that the method dates back to the 
original work of Bauer |^| and Muir [^lj in the 1870's. 

The Bauer-Muir transform of a continued fraction bo{z)+K (a n (z) / b n (z)) with 
respect to a sequence of complex numbers {w n }'^' =0 is defined as the continued 
fraction do(z) + K (c n (z) j d n {z)) , whose "classical" approximants S n (0, z) are equal 
to the modified approximants S n (w ni z) of the original continued fraction. The 
transformed continued fraction exists and can be calculated as 

d = bo + w , Ci=Ai, di=b 1 +wi, (3.37) 
Ci = ai_i^_i : di = bi + Wi - Wi-2qi-i, i > 2, 
Ai = ai — Wi-i(bi +Wi), qi = \i+i/\i i > 1, 

if and only if Aj ^ for i = 1, 2, . . . . 



3.3.2 Pincherle's theorem 

Continued fractions and three-term recurrence relations are intimitaly connected. 
The existence of a minimal solution of a ( 3.15| ) three-term recurrence relation is 



strongly related to the convergence of a continued fraction constructed from the 
coefficients of the recurrence relation. This connection is revealed by Pincherle in 
1894 ||^]. Below we state Pincherle's theorem without its proof which can be found 

in mm. 

Theorem 2 (Pincherle's Theorem). Let be a sequence of complex 

numbers with a n ^ for n= 1,2,... 
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(A) The three-term recurrence relation 

x n+ i = b n x n + a„x„_i, n= 1,2,3,..., (3.38) 

with coefficients a n andb n has a minimal solution {m n } if and only if the continued 
fraction 

b=i V K J h+b 2 + +b n + K ' 1 

converges. 

(B) Provided ( 3.38[ ) Zias a minimal solution {m n }, then for N = 0, 1,2, . . . , 



m jv+i _ j^oo / a«+Af \ _ «i+jv fl2+jv _ a n+jv _ (3 40) 

TOJV ™ _1 \ ^n+Af / + b 2 +N + + b n+N + 

Remarks: 



Equation ( 3.4C ) asserts that the ratio of two arbitrary successive elements of 



the minimal solution can be calculated by a continued fraction. 

• The connection between continued fractions and three-term recurrence rela- 
tions provides the link between continued fractions and special functions (i.e. 
hypergeometric functions and orthogonal polynomials). 

• The genious Indian mathematican Ramanujan left behind, written in his 
notebook [Q, many interesting and original contributions to modern math- 
ematics among which several dealt with the theory of continued fractions. 
Unfortunately he did not give proofs of his ideas. It is worth noting that 
mathematicians have found Pincherle's theorem as a useful tool to prove 
some of Ramanujan's formulae. 



Chapter 4 



Continued fraction representation 
of Green's operators 

In this chapter we present a rather general and easy-to-apply method for calcu- 
lating discrete Hilbert space basis representation of those Green's operators that 
correspond to Hamiltonians having infinite symmetric tridiagonal (i.e. Jacobi) ma- 
trix form. The procedure necessitates the evaluation of the Hamiltonian matrix on 
this basis. The analytically calculated elements of the Jacobi matrix are used to 
construct a continued fraction, which is utilized in the evaluation of the Green's 
matrix. The constructed continued fraction representation of the Green's operator 
is convergent in the bound state energy region and can be continued analytically 
to the whole complex energy plane. 

Our method of calculating Green's matrices ensures a complete analytic and 
readily computable representation of the Green's operator on the whole complex 
plane. Furthermore this is achieved at a very little cost: in practice only the matrix 
elements of the Hamiltonian are required. 

After the expose of the method the continued fraction representation of specific 
Green's operators are given. The Coulomb Green's operator, relativistic Green's 
operators, the Green's operator corresponding to the D-dimensional harmonic os- 
cillator and the Green's operator of the generalized Coulomb potential is consid- 
ered. The convergence and the analytic continuation of the continued fraction is 
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illustrated with the example of the Coulomb Green's operator. The numerical ac- 
curacy of the method is demonstrated via the calculation of the relativistic energy 
spectrum of hydrogen-like atoms. 



4.1 The method 

In order to determine a matrix representation of the Green's operator, first a suit- 
able basis has to be defined. Let us consider the set of states {|i)} and {\i}}, with 
i = 0, 1, 2, . . . , which form a complete biorthonormal basis, i.e. 



OO CO 

i = £l*M = 5»<i| 

and render the E — H operator symmetric tridiagonal 

( ' ■ ^ 

Ji>i = (i'\(E - H)\i) = 



(4.1) 



V 



(4.2) 



/ OO X c 



Here E is a complex parameter and H denotes the Hamilton operator. According 



to Eq. (4.2) we say that the Hamiltonian exhibits a Jacobi matrix structure. 

The Green's operator corresponding to the Hamiltonian H satisfies the operator 
equation 



1 = G(E)(E — H). 



(4.3) 



In the (4.1) discrete biorthonormal basis representation this relation takes the form 
of the following matrix equation 



Sji = £<i|G(^)|?')<i'|(S-ff)|i> Y. (: . 



(4.4) 



i'=0 



Since the basis is chosen such that E — H possesses a tridiagonal matrix represen- 
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tation, the Green's matrix appears as the inverz of a symmetric infinite tridiagonal 
matrix. 



In Section 3.1 the inverz of an infinite symmetric tridiagonal matrix was studied 
and a formula for its n-th leading submatrix was derived. Let us denote the n-th 
leading submatrix of the infinite Green's matrix by G^ . According to Theorem 
|], G^ can be written as 



G 



(n) 



Gon^ 



nn+1 



G()n 



Eq. ([4.5D asserts that the Green's matrix elements can be calculated from the 



(4.5) 



Jacobi matrix elements provided the Go n +i/Go n quotient is at our disposal. 



On the other hand, from the tridiagonality of E — H and from Eq. (4.4) it 
automatically follows that the matrix elements Gji = (j\G(E)\i) satisfy a three- 
term recurrence relation 



Gji-iJi- 



GjiJa 



Gji+iJi+n- 



(4.6) 



Therefore Go„+i/Go n represents a ratio of two consecutive elements of a solution 
of a three-term recurrence relation. It is known, that the solutions of three-term 
recurrence relations span a two-dimensional space and a special type of solution, 
called the minimal solution, is intimately connected to continued fractions, (see 



Section 3.2 and 3.3). According to Pincherle's Theorem, if the Green's matrix 



elements represent a minimal solution of the (4.6) three-term recurrence relation, 
then the ratio of its two consecutive elements is given in terms of a convergent 
continued fraction 



Gon^ 



Gon 



(4.7) 



U n +1 



dn 



U n +2 



+ 1 



d n +2 + 



with coefficients 
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The Green's matrix elements Gij and the ( |4.6| ) three-term recurrence relation, 
by the (4.3) definition, obviously depend on the (complex) energy parameter E. 
First we show that there is a region of the complex S-plane where the physically 
relevant solution of the (4.6) recurrence relation for the Green's matrix is the min- 
imal one, which makes possible a convergent continued fraction representation of 
the Green's operator on this region by utilizing Eqs. (|4.5D and (4.7). Afterwards the 



analytic expression of the convergent continued fraction is extended to other do- 
main of the .E-plane, where the physical Green's matrix is not the minimal solution 
of the recurrence relation. 



On the IRE < region of the complex plane and in case of short-range potentials 
the coordinate space representation of the Green's operator can be constructed as 

@ 

G(r,r',k) = W (fc,r<)// +) (fc,r>)/^(fc), (4.9) 

where ipi(k,r) is the regular solution, f[ (k,r) is the Jost solution, T(k) is the 
Jost function and k is the wave number. The Jost solution is defined by the relation 



lim e Tlfcr / ) '{k,r) = 1. (4.10) 

r— *oo 

Let us define a "new" Green's function as 

G(r, r', k) = r< )Mk, r>)/^(A), (4.11) 

where // is a linear combination of /j and /[ ' . If ^RE < /j is exponentially 
decreasing and f\ ' is exponentially increasing. Thus, for any G we have 

lim ^ r ' ,k) = 0, if IRE < 0. (4.12) 

r'->oa G(r, r', k) 



We note, that both G and G satisfy the defining equation Eq. (|4,3D, but only G of 



Eq. (4.£) is the physical Green's function. The above considerations, with a slight 



modification in Eq. (4.1C), are also valid for the Coulomb case. An interesting result 



of the study of Ref. |54| is that the Green's matrix from Jacobi-matrix Hamiltonian, 
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has an analogous structure to Eq. ( [i~S| ) 

G H ,(k) = ( W ),<(fc)(// +) ) l >(fc)/^(fc), (4-13) 

where (ipi)i{k) = (tpi(k)\l) and (// +) ) 4 (fc) = (// +) (fc)|T). Similarly, we define 
(fiUk) = (fi(k)\i) and 

G w (k) = (w)i<(*)(/0i>(*)/-F(*)- (4.14) 

On the 3?i? < region of the complex plane as r — > oo fi(k,r) exponentially 
domin 
holds 



dominates over fi(k,r), thus for their L 2 representation the following relation 



1™ ^FTVTTT^ = °> if < 0- (4.15) 
This implies a similar relation for the Green's matrices 

l im = 0, if m < 0. (4.16) 

So, we can conclude that in the $IE < region of complex _E-plane the physically 



relevant Green's matrix Gw appears as the minimal solution of the (4J3) recurrence 
relation. Therefore our Green's matrix for bound-state energies can be determined 
by a convergent continued fraction. 

In other regions of the complex energy plane (i.e. the region of scattering en- 
ergies and resonances) the continued fraction fails to converge and the recurrence 
relation does not have a minimal solution. On the other hand, GW+i/GW in (fO) 
is an analytic function of the complex energy parameter, and there is a domain of 
the complex plane where this function is represented by a convergent continued 
fraction, thus values on other domains can be obtained by the analytic continua- 



tion of a continued fraction (see Section 3.3.1). 



Let us suppose we have a limit 1-periodic continued fractions, that is the Ui 
and di coefficients in (|4.8|) have the following limit property 



lim m = u and lim di = d. 

i — >oo i — >oo 



(4.17) 



4.1. THE METHOD 



33 



In this case the continued fraction ( [4,7] ) takes the form 



G 0n u n+1 

d n H 



(4.18) 



d n +i + 



u 

d 



The w tail of this continued fraction satisfies the implicit relation 



w = — U — , (4.19) 

a + w 



which is solved by 



w± 



-d/2± v /(d/2) 2 + u . (4.20) 



Replacing the tail of the continued fraction by its explicit analytic form w± , we 
can speed up the convergence and, which is more important, we can perform an 



analytic continuation. This implies the usage of w± — w±(z) in the (3.25) formula 
for the S n (w, z) continued fraction approximants also at those z values, where the 
continued fraction fails to converge. 

The w+ choice gives an analytic continuation to the physical sheet, while u;_, 
which also converges, gives an analytic continuation to the unphysical sheet. This 
observation can be taken by considering the formula for the Green's operator on 
the physical sheet pGp 

(1\G(E + iO)|T) - {i\G{E - iO)|T) = -27ri(7|V(£0)<^(£)|T), (4.21) 



where ip(E) is the scattering wave function, and utilizing the ( |2.1l| ) analytic prop- 
erties of Green's operators. Then we can readily obtain that the imaginary part of 
(i\G(E + i0)\i) should be negative and this condition can only be fulfilled with the 
choice of w+ . 



We have shown that Eqs. (4.5) and ( |4.7D together with the theory of analytic 



continuation of a continued fraction supply a complete discrete Hilbert space repre- 
sentation of the Green's operator on the whole complex energy plane. The Green's 
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matrix can be obtained from the Jacobi Hamiltonian for arbitrary complex energies 
by simply evaluating a continued fraction and performing a matrix inversion. 



4.2 D-dimensional Coulomb Green's operator 

Here we use the Coulomb-Sturmian basis and show that on this basis the D- 
dimensional Coulomb Hamiltonian possesses a Jacobi-matrix structure. Then the 
analytically derived J-matrix elements are utilized, following the method of the 
previous section, to calculate the continued fraction representation of the Green's 
operator. The convergence of the continued fraction and the technique of the 
analytic continuation is demonstrated in practice. 

Let us consider the .D-dimensional radial Coulomb Hamiltonian 

etc h 2 f d 2 1/ D-3U, D-l\\ Ze 2 



2m \dr 2 r 2 \ 2 J \ 2 J J r 

where m, I, e and Z stand for the mass, angular momentum, electron charge and 
charge number, respectively. (See, e.g. Ref. JHjj] and references.) The bound state 
energy spectrum is given by 

mZ 2 e 4 

En > 4 = " 2h 2 (n r + l+°f±) 2 ' (4 - 23) 
and the corresponding wave functions are 

(4.24) 

where we used the notation ao = ((n r + I + : ^r-)ro)~ 1 and tq = h 2 /(2mZe 2 ). 

The Coulomb-Sturmian (CS) functions, the Sturm-Liouville solutions of the 
Hamiltonian ( 4.22 ), appear as 

1 /2 

«M 6 ' r >= (fj^Flj) exp(-6r)(2&r)' + ^4 2i+D - 2) (26r), (4.25) 
where b is a scale parameter, n is the radial quantum number and L„ denotes 
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the generalized Laguerre polynomials |Q. The CS functions of ( 4.25 ) are the gen- 
eralizations of the corresponding CS functions for the three-dimensional case ||57|| . 
Introducing the notation (r|n) = <j> n i(b,r) and (r\n) = <f) n i(b,r)/r for the CS func- 
tion and its biorthonormal partner respectively, we can express the orthogonality 
and completeness of these functions as 



(n\ri) = 5 n 

OO OO 

l = X;|n)<n| = 53|n)(f 



(4.26) 
(4.27) 



confirming that they form a discrete biorthonormal basis in the sense of (4.1). The 
overlap of two CS functions can be written in terms of a three-term expression 



(n\n') = (2b)- 1 + <W (2n + 21 + D - 1) 

- S nn '-i((n + l)(n + 21 + D — 1)) 1/2 
-5 nnl+ i(n{n + 2l + D-2)) 1 ' 2 . 

A similar expression holds for the matrix elements of H 



(n\H c \n') 



4m 



/ 2 

+ 5 nn , [2n + 2l + D -I - 

V r ob 

+ <W-i((n + l)(n + 21 + D - l)) 1 / 2 
+ 5 nn , +1 (n(n + 2l + D-2)) 1 / 2 . 



(4.28) 



(4.29) 



Let us denote the Coulomb Green's operator as G (E) = (E — i/ c ) _1 and 
determine its CS matrix elements Gfj = (i\Gf\j) by applying the general method 
described previously. The starting point in this procedure is the observation that 
the Jij = (i\(E — H c )\j) matrix possesses an infinite symmetric tridiagonal i.e. 
Jacobi structure, and the nonzero elements of this tridiagonal matrix are obtained 
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immediately from Eqs. ( 4.28 ) and ( 4.29 ) as 



h 2 

J u = (2i + 2l + D -l)(k 2 - b 2 ) --Ze 2 (4.30) 

Amb 

Ja-i = -[i(i + 2l + D-2)] 1 / 2 (k 2 + b 2 )^, 

where k = (2mE/h 2 ) 1/2 is the wave number. 

Then the n-th leading submatrix G^™" 1 of the infinite Green's matrix is repre- 
sented by 

Gy (n) = [Jij + 5 jn 5 in Jnn+1 F}' 1 , (4.31) 

where F is a continued fraction 



F = , (4.32) 

U n +1 

On H 

U n +2 

d n +l 



d n +2 H 

with coefficients 

Ui = —Ji,i-i/ Ji,i+l, di = —Ji,i/ Ji,i+1 ■ (4.33) 

The above continued fraction F, as it stands, is only convergent for negative 
energies, but since it is a limit-1 periodic continued fraction, i.e. its coefficients m 
and di possess the limit properties 

u = lim Ui = — 1 (4.34) 

i — >oo 

d = lim di = 2(k 2 - b 2 )/(k 2 +b 2 ) , 

i — >oo 

it can be continued analytically to the whole complex energy plane by replacing its 
tail with 



w± = (b±ik) 2 /(b 2 + k 2 ) , 



(4.35) 
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according to Section 4.1. This way Eq. ( 4.31 ) provides the CS basis representation 
of the Coulomb Green's operator on the whole complex energy plane. 

We note here that the Coulomb-Sturmian representation of the Coulomb Green's 
operator Gfj has already been calculated |g5|, ^?], ^8|, by evaluating a complicated 



integral of the coordinate space Green's operator and the ( 4.25 ) CS functions. This 
integral, for example in the case of Gg = (0|G C |0) could be performed analytically 
leading to the result 



Amb 



K 2 (b-ik) 2 l + (D-l)/2 + ij 



x 2 Fi -I - 



D-3 , D 

— + vy,l;l + — 



•17' 




(4.36) 



containing the 2-F1 hypergeometric function |Q. Afterwards Gq together with a 
three-term recurrence relation could be applied in order to calculate other matrix 
elements recursively. 

We point out that with the choice of Z — the D-dimensional Coulomb Hamil- 



tonian (|4.22| ) reduces to the D-dimensional kinetic energy operator and our formulas 
provide the CS basis representation of the Green's operator of the free particle as 
well. 



4.2.1 Convergence of the continued fraction 



Below we demonstrate the convergence and the numerical accuracy of the ( 4.31 ) 
continued fraction representation of the Coulomb Green's operator. We calculate 
the Gqq(E) matrix element of the D-dimensional Coulomb Green's operator for 
/ = and D = 3 case at E = (—100,0) and E = (1000,1) energy values from 
the bound and scattering state regions of the complex _E-plane. For comparison, 



the numerical value of Gq (E) calculated by the analytic expression (4.36) is also 
quoted here and denoted by Ga- 



The continued fraction ( 4.32 ) has been evaluated by calculating its n-th ap- 
proximants with respect to different {w n } series. For bound state energies the 
convergence of the continued fraction with respect to the different choice of {w n }, 
while for the scattering case the effect of the analytic continuation and the Bauer- 
Muir transformation has been examined. 
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In Table 4.1 we can observe excellent convergence of the continued fraction 
to the exact value in all cases. In case of $IE < the Green's matrix elements 
represent the minimal solution of a three-term recurrence relation, thus due to 
Pincherle's theorem the continued fraction is convergent. The different choice of 
{w n }, e.g. we took the w n — 0, w n = w+ and w n — w- choices, influences only the 
speed of convergence. 



In Table 4.2 the continued fraction approximants of Gqq(E) are shown for a 
scattering energy case. In the region of 5tE > 0, in complete accordance with 
Pincherle's theorem, the original continued fraction (4.32) diverges, only its analytic 
continuations are convergent. We recall here that w n — w+ provides analytic 
continuation to the physical, while w n — W- to the unphysical sheet. In Table 
4.2 only approximants with respect to w+ are given. However, as the first column 
shows, the convergence is rather poor. This can be considerably improved by the 
repeated application of the Bauer-Muir transformation (see Section 3.3.1 ). In fact, 
as can be seen in the last column, an accuracy similar to the bound state case can 
be easily reached here with e.g. an eightfold Bauer-Muir transform. 



Examining Table 4.1 and 4.2 we can draw the conclusion that the general and 
easily computable continued fraction method of Section 11 provides a convergent, 
numerically stable and accurate representation on the whole complex plane. 



4.2.2 Numerical test 

We can immediately test the analytic properties of G c by determining the ( [4.23] ) 
eigenvalues of the attractive Coulomb interaction in three-dimension as the poles of 



the Green's matrix (see Section g). Figure 4T shows &et[G c (E)}, the determinant 
of the Green's matrix as the function of the E energy parameter. The poles coincide 
with the exact Coulomb energy levels up to machine accuracy. We stress that, from 
the point of view of determining the energy eigenvalues, the rank of the matrix and 
the specific choice of the CS basis parameter are irrelevant. An arbitrary-rank 
matrix representation of the Coulomb Green's operator exhibits all the properties 
of the system and our Green's matrix contains all the infinitely many eigenvalues. 
This is especially interesting if we compare with the usual procedure of calculating 
eigenvalues of a finite Hamiltonian matrix of rank N, which could only provide 
an upper limit for the N lowest eigenvalues. Our procedure does not truncate 
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n 


G (o) 


E = (-100,0) 

G (w + ) 




1 


(-5.44922314793965,0) 


(-5.59142801316938,0) 


[-0.92906408986331,0) 


2 


(-5.54075476366523,0) 


(-5.56131039101044,0) 


[-4.70080363351349,0) 


3 


(-5.55501552420656,0) 


(-5.55812941271530,0) 


[-5.39340957492282,0) 


4 


(-5.55726787304507,0) 


(-5.55775017508704,0) 


[-5.52662100417805,0) 


5 


(-5.55762610832912,0) 


(-5.55770176067796,0) 


[-5.55192403535264,0) 


6 


(-5.55768333962797,0) 


(-5.55769530213874,0) 


[-5.55663951922343,0) 


7 


(-5.55769251168083,0) 


(-5.55769441374319,0) 


[-5.55750389878413,0) 


8 


(-5.55769398510141,0) 


(-5.55769428874846,0) 


[-5.55766025755765,0) 


9 


(-5.55769422223276,0) 


(-5.55769427085427,0) 


[-5.55768824220786,0) 


10 


(-5.55769426045319,0) 


(-5.55769426825710,0) 


[-5.55769320762502,0) 


11 


(-5.55769426662100,0) 


(-5.55769426787592,0) 


[-5.55769408236034,0) 


12 


(-5.55769426761735,0) 


(-5.55769426781946,0) 


[-5.55769423553196,0) 


13 


(-5.55769426777843,0) 


(-5.55769426781103,0) 


[-5.55769426221577,0) 


14 


(-5.55769426780450,0) 


(-5.55769426780976,0) 


[-5.55769426684377,0) 


15 


(-5.55769426780872,0) 


(-5.55769426780957,0) 


[-5.55769426764335,0) 


16 


(-5.55769426780940,0) 


(-5.55769426780954,0) 


[-5.55769426778103,0) 


17 


(-5.55769426780951,0) 


(-5.55769426780954,0) 


[-5.55769426780465,0) 


18 


(-5.55769426780954,0) 




[-5.55769426780870,0) 


19 


(-5.55769426780954,0) 




[-5.55769426780939,0) 


20 






[-5.55769426780950,0) 


21 






[-5.55769426780954,0) 


22 






[-5.55769426780954,0) 




Ga 


= (-5.55769426780954, 





Table 4.1: Convergence of the continued fraction for the Goo element of the Green's 
matrix at $IE < 0. The first, second and third column contain approximants of 
the continued fraction with w n — 0, w n = w + and w n — respectively. For 
comparison we also give the Ga ( 4.36 ) analytic result. All the Goo values are 
scaled with 10 2 . 
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n 


G^+^O) 




G (»+)( 1 ) 


E = (1000, 1) 


G (w +\8) 


1 


l J..U i U 5 


-U.U IOJ 


(1 


8079 


-u.ozyo j 


C n ii 9Q^4zl 




(-0.2321154,-0.073120618) 


5 


(1 D74 




(1 


1225 






1 R3494Q31 *l 


(-1.4408861,-0.350899497) 


10 


(1.198, 


-0.325) 


(1 


1425, 


-0.3162) 


(1.13445656, 


-0.32244006) 


(1.20667672,0.237375310) 


15 


(1.110, 


-0.346) 


(1 


1497, 


-0.3353) 


(1.14598003, 


-0.33019962) 


(1.14702383,-0.332562329) 


20 


(1.160, 


-0.307) 


(1 


1415, 


-0.3287) 


(1.14512823, 


-0.33023791) 


(1.14511731,-0.330140243) 


25 


(1.141, 


-0.354) 


(1 


1478, 


-0.3298) 


(1.14523860, 


-0.33015825) 


(1.14523597,-0.330179581) 


30 


(1.139, 


-0.312) 


(1 


1437, 


-0.3311) 


(1.14522511, 


-0.33018993) 


(1.14522395,-0.330182552) 


35 


(1.157, 


-0.341) 


(1 


1458, 


-0.3290) 


(1.14522377, 


-0.33017859) 


(1.14522539,-0.330181124) 


40 


(1.131, 


-0.325) 


(1 


1451, 


-0.3312) 


(1.14522642, 


-0.33018226) 


(1.14522527,-0.330181563) 


45 


(1.158, 


-0.327) 


(1 


1448, 


-0.3294) 


(1.14522458, 


-0.33018138) 


(1.14522524,-0.330181440) 


50 


(1.135, 


-0.337) 


(1 


1457, 


-0.3305) 


(1.14522559, 


-0.33018133) 


(1.14522527,-0.330181470) 


55 


(1.149, 


-0.320) 


(1 


1446, 


-0.3300) 


(1.14522512, 


-0.33018161) 


(1.14522525,-0.330181465) 


60 


(1.145, 


-0.340) 


(1 


1456, 


-0.3300) 


(1.14522528, 


-0.33018135) 


(1.14522526,-0.330181464) 


65 


(1.140, 


-0.322) 


(1 


1449, 


-0.3304) 


(1.14522527, 


-0.33018153) 


(1.14522525,-0.330181466) 


70 


(1.152, 


-0.334) 


(1 


1453, 


-0.3298) 


(1.14522523, 


-0.33018143) 


(1.14522526,-0.330181465) 


75 


(1.137, 


-0.329) 


(1 


1452, 


-0.3304) 


(1.14522528, 


-0.33018147) 


(1.14522526,-0.330181466) 


80 


(1.152, 


-0.327) 


(1 


1450, 


-0.3296) 


(1.14522524, 


-0.33018146) 


(1.14522526,-0.330181465) 


85 


(1.140, 


-0.335) 


(1 


1454, 


-0.3302) 


(1.14522527, 


-0.33018145) 


(1.14522526,-0.330181465) 


90 


(1.147, 


-0.323) 


(1 


1450, 


-0.3301) 


(1.14522525, 


-0.33018147) 


(1.14522526,-0.330181465) 


95 


(1.146, 


-0.336) 


(1 


1453, 


-0.3300) 


(1.14522526. 


-0.33018145) 


(1.14522526,-0.330181465) 












Ga = 


(1.14522526, - 


-0.330181465) 





Table 4.2: Convergence of the continued fraction for the Goo element of the Green's 
matrix at ?R.E > 0. All columns contain approximants of the analytic continuation 
of the continued fraction with respect to w+. While the first, second, third and 
fourth column contain approximants without, with one-fold, with five-fold and with 
eight- fold Bauer-Muir transform, respectively. For comparison we also give the Ga 
( 4.36 ) analytic result. All the Goo values are scaled with 10 2 . 
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the Coulomb Hamiltonian, since all the higher Jij matrix elements are implicitly 
contained in the continued fraction. 

In order to have a more stringent test we have performed the contour integral 

2^£ dEGoo(E). (4.37) 

If the domain surrounded by C does not contain any pole, then 1(C) = 0. If this 
domain contains a single bound state pole, 1(C) = (0|^>)(^|0) must hold, while if C 
circumvents the whole spectrum then 1(C) = (0|0) is expected. With appropriate 
selection of Gauss integration points we could reach 12 digits accuracy in all cases. 
This indicates that the calculation of the Green's matrix from a J-matrix via the 
continued fraction method is accurate on the whole complex plane. 

4.3 Relativistic Coulomb Green's operators 

In this section we specify our method for relativistic Coulomb Green's operators: 
the Coulomb Green's operator of the Klein-Gordon and of the second order Dirac 
equations. The latter is physically equivalent to the conventional Dirac equation 
but seems to have several advantages from the mathematical point of view. For 
details see Ref. j58| and references therein. 

The Hamiltonian of the radial Coulomb Klein-Gordon and second order Dirac 
equations are shown to possess an infinite symmetric tridiagonal matrix structure 
on the relativistic Coulomb-Sturmian basis. This allows us to give an analytic 
representation of the corresponding Coulomb Green's operators in terms of con- 
tinued fractions. The poles of the Green's matrix reproduce the exact relativistic 
hydrogen spectrum. 

It is noted here that the Coulomb-Sturmian matrix elements of the second 
order Dirac equation has already been obtained by Hostler |58| via evaluating 
complicated contour integrals. Our derivation, however is much simpler, it relies 
only on the Jacobi-matrix structure of the Hamiltonian, and the result obtained is 
also better suited for numerical calculations. In Hostler's paper the results appear 
in terms of T and hypergeometric functions, while our procedure results in an easily 
computable and analytically continuable continued fraction. 

The radial Klein-Gordon and second order Dirac equation for Coulomb inter- 
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1 


1 


1 


1 






1 


1 


1 


1 







-0.2 -0.15 -0.1 -0.05 
E (atomic units) 



Figure 4.1: The determinant of a 3 x 3 three-dimensional Coulomb Green's matrix 
det[G c (E)] as the function of the energy E for I = 1. The bound states of the 
Coulomb problem are located at energies where the vertical lines cross the hori- 
zontal axis. (These lines are shown only for demonstrative purposes, they do not 
correspond to the values of det[G (E)].) Atomic units ofm = c = ft=l and 
Z = — 1 were used. For the sake of clarity only the first 6 energy levels are shown. 
These are located at E n = — l/[2(n + I + l)] 2 with n < 5 according to (4.23). 
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action are given by 



h u ID = o, 



(4.38) 



where 



2aZ E d 2 u(u + l) 
he r dr 2 r 2 



(4.39) 



Here /U = mc/h, a — e 2 /hc, m is the mass and Z denotes the charge. For the 
Klein-Gordon case u is given by 



u = -\ + J\+ 1(1 + 1) -(Za) 2 , 



(4.40) 



and in the case of the second order Dirac equation for the different spin states we 
have 



(4.41) 



The relativistic Coulomb Green's operator is defined as the inverse of the Hamil- 
tonian H u 



H U G U G U H U l u , 
where l u denotes the unit operator of the radial Hilbert space H. u . 



(4.42) 



In complete analogy with the non-relativistic case we can define the relativistic 
Coulomb-Sturmian functions as solutions of the Sturm-Liouville problem 



.* + ^ + S<S±li-2Sfit±5±a")«i,(r)-0, ,4-4,, 



where 77 is a real parameter and n = 0, 1, 2, . . . ,00 is the radial quantum number. 
They take the form 



(r \nu;r)} = S%. v (r) 



IV. 



(n + 2u+l)! 



(2r,r) u+1 e-^Ll u+1 (2r,r), (4.44) 
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where L is a Laguerre-polinom. The relativistic Coulomb-Sturmian functions, 
together with their biorthonormal partner (r \nu; rj) — l/r ■ (r \nu; n), form a basis: 
i.e., they are orthogonal 

(mi; rj \mu; r/) — (nu; rj \mu; rf) = 8 nm , (4.45) 

and form a complete set in Ti. u 

\nu; rj) (nu; j]\ — \nu;rj) (nu;r)\ = l u . (4.46) 

n=0 n=0 

A straightforward calculation yields 



1 



(nu;rj\mu; rj)=— 5 nm (2u + 2n + 2) - 5 n m-is/ (n + l)(n + 2u + 2) 
Zrj L 



-5 nm+ i y/n(2u + n + l) 



(4.47) 



Utilizing this relation and considering Eq. ( 4.43| ) we can easily calculate the Coulomb- 
Sturmian matrix elements of H u , 

H nm ■= (nU] T]\ H u \mU\ 1]) = 

. (2azE . n , .(E/hc) 2 - Li 2 + n 2 

+ 5 nm [— 2(u + n+l)n + 2(u + n+l) y 



he 2r\ 

fiElhc) 2 - f i 2 +n 2 n - r A (4-48) 

- <W-1 f 1 ' ' ^ '-^{n + l){n + 2u + 2)\ 

f(E/hc) 2 - u 2 +n 2 i—. A 

- S nm+1 i K 1 ' ; '-^n(n + 2u+l) \ , 

which happen to possess a Jacobi-matrix structure. 



Now the Coulomb-Sturmian matrix elements of the relativistic Coulomb Green's 
operators 

( G «)„m = (««; ri\G u \mu; n) (4.49) 



corresponding to Hamiltonians (4.39), can be straightforwardly determined by us- 
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ing the continued fraction method of Equations (Ejj) and ( [4,7] ) with the coefficients 

Ha- 1 , Hi,, 



H, 



di = — 



ii+l 



H 



ii+l 



(4.50) 



Again, the continued fraction representation convergent for bound-state energies 
and can be continued analytically to the whole complex energy plane. 



4.3.1 Relativistic energy spectrum 



In Table 4.3 we demonstrate the numerical precision of our Green's matrix by 
evaluating the ground and some highly excited sates of relativistic hydrogen-like 
atoms, which, in fact, correspond to the poles of the Dirac Coulomb Green's matrix. 
In particular, the poles of the determinant of ( 4.49|) were located. Here we repeat 
again that irrespective of the rank N of the Green's matrix the poles should provide 



the exact Dirac results. In Table |13| we have taken 2x2 matrices. Indeed, the 
results of this method, E c {, agree with the exact one in all cases, practically up to 
machine accuracy, this way making possible the study of the fine structure splitting. 





energy levels 


E ( , 






hydrogen 


lSl/2 


-0.5000066521 


-0.5000066521 


-0.5 


2Pl/2 


-0.1250020801 


-0.1250020801 


-0.125 


2P 3/2 


-0.1250004160 


-0.1250004160 


-0.125 


50P 1/2 


-0.0002000002 


-0.0002000002 


-0.0002 


50P 3/2 


-0.0002000001 


-0.0002000001 


-0.0002 


uranium 


lSl/2 


-4861.1483347 


-4861.1483347 


-4232 


100D 3/2 


-0.4241695002 


-0.4241695002 


-0.4232 


100D 5/2 


-0.4238303306 


-0.4238303306 


-0.4232 



Table 4.3: Energy levels of hydrogen-like atoms in atomic units. E c ^ is the relativis- 
tic spectrum calculated via continued fraction, E pj and -Eg are textbook values of 
the relativistic Dirac and the non-relativistic Schrodinger spectrum, respectively. 



4.4 D-dimensional harmonic oscillator 

The Hamiltonian of the D-dimensional harmonic oscillator reveals a Jacobi matrix 
representation in the harmonic oscillator basis. Therefore the analytically calcu- 
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lated Jacobi-matrix elements can be utilized as the input of the continued fraction 
method for calculating harmonic oscillator basis representation of the Green's op- 
erator corresponding to the harmonic oscillator potential. 



The radial Hamiltonian describing the D-dimensional harmonic oscillator prob- 
lem has the form 



H HO = 



h 2 ( d 2 1 /, £>-3\ /, D-l\\ 1 



1 + 



+ 2 muj r 



where u> is the harmonic oscillator parameter. The energy eigenvalues are 

E n i = hu [ 2n + I + D 



(4.51) 



(4.52) 



and the corresponding eigenfunctions can be written as 

1/2 



(r|oj, nl) — v 4 



i / 2T(n + 1) 



T(n + l + §) 



e X p(_| r 2)( ur 2 ) i+^ i ^+*- 1 ) (w 2 )j (4 _ 53) 



where v = mu>/H. The harmonic oscillator functions ( 4.53 ) with fixed u are or- 
thonormal and form a complete set in the usual sense 



(u; , n 1 l\ui , nl) = S mv 

oo 

1 = \^,nl){w,nl\ 



(4.54) 
(4.55) 



The harmonic oscillator Hamiltonian (4.51) with parameter u> on the basis of the 
harmonic oscillator functions with different parameter to' takes a Jacobi-matrix 
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form 



(uj',nl\H HU (uj)\uj',n'l) 



+ S nn ' h 



2uj' 



2ri + l + Y 



2 /2 

UI — LO 



2uj' 



n' \n' + 1 + 



D 



1/2 



Snn'+iti 2uj , ( (n' + 1) ( n 1 + I + — 



1/2 



(4.56) 



The general method of Section [4^1] requires the knowledge of the matrix elements 
Jij = (uj',U\E - H ho (uj)\lj' ,jl) which readily follows from (|4,54[) and ( |4.56| ) ac- 
cording to 



Jij = Ex Sij - (uj',il\H HO (uj)\uj',jl) 



(4.57) 



The calculation of the Green's matrix Gij(E) = (a/, iZ|G(£J, ji) goes similarly 



to the previous sections making use of formulae ( |4.5| ) and (4.7). 

It is impossible to overestimate the importance of the harmonic oscillator in 
theoretical physics. Here I would like only to mention one exotic topic, the physics 
of anyons, which are quantum mechanically indistinguishable particles following 
fractional statistics, where the harmonic oscillator potential plays an important 
role [El. 



4.5 The generalized Coulomb potential 

Quantum mechanical models and practical calculations often rely on some exactly 
solvable models like the Coulomb and the harmonic oscillator potentials. The 
actual example we consider here is the generalized Coulomb potential introduced by 
Levai and Williams |50| , which is the member of the Natanzon confluent potential 
class J6l]]. This potential is Coulomb- like asymptotically, while its short-range 
behavior depends on the parameters: it can be finite or singular as well at the origin. 
Its shape therefore can approximate various realistic problems, such as nuclear 
potentials with relatively flat central part, or atomic potentials that incorporate 
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the effect of inner closed shells by a phenomenological repulsive core. Another 
interesting feature of the D-dimensional generalized Coulomb potential is that it 
contains the Coulomb and harmonic oscillator potentials as limiting cases, this way 
providing a smooth transition between the Coulomb and the harmonic oscillator 
problems in various dimensions. 

More and more interactions can be modelled by making advantage of the rather 
flexible potential shapes offered by exactly solvable potentials. Virtually all quan- 
tum mechanical methods rely in some respect on analytically solvable potentials. 
Very often their wave function solutions are used as Hilbert space bases. More 
powerful methods can be constructed if we select a basis which allows the exact 
analytical calculation of the Green's operator of an analytically solvable potential. 

In this section we show that an appropriate Sturm-Liuville basis can be defined 
on which the matrix elements of the Hamiltonian exhibit a Jacobi matrix. The 



corresponding Green's matrix then follows from the method of Section 4.1 



4.5.1 The potential 



Let us consider the radial Schrodinger equation in D spatial dimensions with a 
potential V(r) that depends only on the radial variable r 



D-l, 



) + v(r) ip(r) = eip(r) , (4.58) 



where v(r) = 2mh 2 V(r) and e = 2mh 2 E. We define the generalized Coulomb 
potential [M in D-dimension as 



«(r)=-^(i + -2-)(i 



D- I s 1 3 C 

2 ) + (/ 2> {l 2Uh{r)(h(r)+6) 



q 



3C 



5C9 



(4.59) 



h(r)+9 16{h(r)+6) 2 16(/i(r) + Of ' 

where h(r) is defined in terms of its inverse function 

i \ 

h 



r = r(h) = C-i 



fltanh" 1 



h + 9 



+ (h(h + 9)y- 



(4.60) 
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The h(r) function maps the [0, oo) half axis onto itself and can be approximated 
with h(r) ~ C^r and h(r) ~ Cr 2 / (49) in the r — ► oo and r — * limits, respectively. 

Bound states are located at 



where 



C 2 



q6 
C 



(4.61) 



(4.62) 



and the bound-state wave functions can be written in terms of associated Laguerre- 
polynomials as 



ipn{r) = C*p„ 



r(n + l) 



T(n + P)(2n + (3 + p n t 



1/2 



x (h(r) +e)i(h(r)) M ^eM-^Kr))L^- 1 \p n h(r)) 



(4.63) 



Potential (4.59) clearly carries angular momentum dependence: its first term 
merely compensates the centrifugal term arising from the kinetic term of the Hamil- 
tonian. Its second term also has r~ 2 -like singularity (due to /i _1 (r)) and it cancels 
the angular momentum dependent term in the two important limiting cases that 
recover the -D-dimensional Coulomb and the harmonic oscillator potentials. The 



third term of (|4.59j ) represents an asymptotically Coulomb-like interaction, while 
the remaining two terms behave like r~ 2 and r -3 for large values of r. The long- 



range behavior of potential ( 4.5S ) suggests its use in problems associated with the 
electrostatic field of some charge distribution. The deviation from the Coulomb 
potential close to the origin can be viewed as if the point-like charge could be 
replaced with an extended charged object. The relevant charge density is readily 
obtained from the potential using 



p{r) 



87rme 



Av(r) 



(4.64) 



In Figures 4.2 and 4.3 we present examples for the actual shape of potential 
( 4.59] ) and the corresponding charge distribution (4.64) for various values of the 
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Figure 4.2: The generalized Coulomb potential for g=0.5, 1.25, 2.5; 0=0.01, 0.1, 1; 
C=l and (3 = 3/2. I — and D = 3 is also implied. In each panel the largest q 
corresponds to the lowest curve. Note the different scales of the horizontal (r) and 
the vertical (v(r)) axes. 



parameters. It can be seen that this potential is suitable for describing the Coulomb 
field of extended objects. It is a general feature of potential ( 4.59 ) that for small 
values of 9 a (finite) positive peak appears near the origin, which also manifests itself 
in a repulsive "soft core", corresponding to a region with positive charge density 
(see Fig. O). 



The Coulomb and harmonic oscillator limits 



The special limits of the generalized Coulomb potential can be realized by specific 



choices of the parameters in Eq. fl4.60| ) : 

The ZJ-dimensional Coulomb limit follows from the 9^0 limit and it is recov- 
ered from Eq. ( fl.59| ) by the = 21 + D - 1 and C~?q = 2mZe 2 /h 2 , choices: the 



third term of (4.59) becomes the Coulomb term, the fifth one vanishes, while the 



other three terms becoming proportional with r 2 cancel out completely. 

In order to reach the oscillator limit one has to take 9 — > oo keeping C/9 = 
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-0.1 










-0.15 










-0.2 










-0.25 











Figu re 4 ,3: The charge distributions corresponding to the potentials displayed in 
Fig. for 9=1 and g=0.5, 1.25, 2.5. (The largest q corresponds to the lowest 
curve.) 



C constant together with the redefinition of the potential ( 4.59Q and the energ 



eigenvalues by adding q/9 to both. This choice simply represents resetting the 
energy scale: e = corresponds to v(r — > oo) for the Coulomb problem, and to 
v(r = 0) for the harmonic oscillator. (Note that the energy eigenvalues also have 
different signs in the two cases.) Besides C/9 = C, the q = q/9 2 parameter also 
has to remain constant in the 9 — > oo transition here. The potential thus adapted 
to the harmonic oscillator limit reads 

v{r) = v(r) + q9 = 



2 A 2 y V 2/ V 2; 4h(r)(l + hip.) 
qh(r) 3C* 1 5C* 1 



1 + 169( 1 + !lirl y 166 ( Mr) 



(4.65) 



The harmonic oscillator potential is recovered from ( 4.65| ) by the (3 = I + D/2 
and Cq = (2muj/h) 2 choice. The two last terms in (4.65) vanish, the first and the 
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second cancel out, while the third one reproduces the harmonic oscillator potential. 
The new form of the energy eigenvalues is 



e„ + q/9 = C{2n + 0) 



2 ' 



i 



2' 



(4.66) 



which indeed, reduces to the e n = (2mw /h)(2n + l + D /2) oscillator spectrum in the 
8 — > oo limit. The wave functions ( 4.63| ) are unchanged, except for the redefinition 
of the parameters. 



4.5.2 The matrix elements of the Green's operator 



We define the generalized Coulomb-Sturmian basis as the solution of the general- 
ized Sturm-Liouville equation. The Sturm-Liouville equation, which depends on n 



as a parameter and corresponds to the generalized Coulomb potential (|4.59|) , reads 



3C 



5ce 



c(/3 -!)(/?-!) 



dr 2 W(h(r) + 6) 2 16(h{r) + 9) 3 4h(r)(h(r) + 6) 
- 2a (3 S \ C C 



-\—+P( n +2 ) Jh(r) 



1 P 



tf>(p,r)=0, (4.67) 



and is solved by the generalized Coulomb-Sturmian (GCS) functions 



(r\n) = 4> n (p,r) 



r(n + l) 
r(n + 0) 



1/2 



{ph(r) + p9) 2 (ph(r)) ^ 1 exp( 



-h(r))L^(ph(r)) . (4.68) 



Here p is a parameter characterizing the generalized Coulomb-Sturmian basis. The 
GCS functions, being solutions of a Sturm-Liouville problem, have the property of 
being orthonormal with respect to the weight function C*(h(r)+0)~ x . Introducing 
the notation (r\n) = <j) n (p,r)C^ (h(r) + O)^ 1 the orthogonality and completeness 
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relation of the GCS functions can be expressed as 

(n'\n) = S n > n 
l = f>>(n|=f>><? 



(4.69) 



n=0 



Analytic calculations yield that both the overlap of two GCS functions and the 
H n / n = (n'\H\n) Hamiltonian matrix possesses a tridiagonal form, therefore the 
matrix elements of the e — H operator also have this feature 



(n\e-H\n') 



e C 2 o 



<W+i (n(n + /3-l))s | 
<W-i((n + l)(n + /?)P 



4q_ 

Cp 



(2n + P) 



Cip 
/ e 



c -¥) 

Cip" 



(4.70) 



Cip 



This means, that similarly to the previous sections the matrix elements of the 
Green's operator in the GCS basis, Gy = (n\G\n), can be determined by using 
continued fractions, as described in Section JO] utilizing the analytically known 
Jacobi-matrix elements of (4.70). 



Chapter 5 

Applications 



The continued fraction method for calculating Green's matrices on the whole com- 
plex energy plane together with methods for solving integral equations in discrete 
Hilbert space basis representation provide a rather general and easy-to-apply quan- 
tum mechanical approximation scheme. 

In the first part of this chapter the continued fraction representation of the 
Coulomb-Sturmian space Coulomb Green's operator (Section 4^) is used for solv- 
ing the two-body Lippmann-Schwinger equation with a potential modelling the 
interaction of two a particles in order to find bound, resonance and scattering 
solutions. 

In the second part of this chapter our Green's operator is applied for solving 
the Coulomb three-body bound state problem in the Faddeev-Merkuriev integral 
equation approach. In particular, the binding energy of the Helium atom is de- 
termined by solving the Faddeev-Merkuriev equations in the Coulomb-Sturmian 
space representation. 



Both solution schemes have been devised by Papp in Refs. [|25|, |27], |28|, 29 1 
and Refs. JL6| , |30[ , respectively. We demonstrate here that the continued fraction 
representation of the Coulomb Green's operator in practice is as good as the original 
one given by Papp in terms of hypergeometric functions. 

The two examples of this chapter are intended to show the importance of the an- 
alytic representation of the Green's operators through the efficiency of the discrete 
Hilbert space expansion method for solving fundamental integral equations. 
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5.1 Model nuclear potential calculation 

In this section we apply the method of Refs. |p5| , p7| , |2S[ , |29|| together with the 
continued fraction representation of the Coulomb Green's operator in order to cal- 
culate bound, resonant and scattering state solutions of a potential problem in a 
unified manner. The particular example we consider here is a potential modelling 
the interaction of two a particles. This example is thoroughly discussed in the ped- 
agogical work |6^| in the context of a conventional approach based on the numerical 
solution of the Schrodinger equation. 

The interaction of two a particles can be approximated by the potential 

V a - a {r) = -Aexp(-/3r 2 ) + — erf( 7 r), (5.1) 

r 

where erf(z) is the error function |56). This potential is a composition of a bell- 
shaped deep, attractive nuclear potential, and a repulsive electrostatic field between 
two extended charged objects. The units used in the Hamiltonian of this system 
are suited to nuclear physical applications, i.e. the energy and length scale are 
measured in MeV and fm, respectively. In these units H/(2m) = 10.375 MeV fm 2 
(m is the reduced mass of two a particles) and e 2 = 1.44 MeV fm. The other 
parameters are A =122.694 MeV, (3 = 0.22 fm" 2 , 7 = 0.75 fm" 1 and Z = 2 (the 
charge number of the a particles) . 



Our radial Hamiltonian Hi containing the model potential (|5.l| ) can be split 
into two terms 

Hi = Hp + Vi. (5.2) 
Here Vi is the asymptotically irrelevant short-range potential and Hp denotes the 



asymptotically relevant radial Coulomb Hamiltonian (4.22). Since the a — a po- 



tential possesses a Coulomb tail, the short-range potential is defined by 

Z 2 e 2 Z 2 e 2 
Vi(r) = V a - a {r) = -Aexp(-/3r 2 ) erfc( 7 r), (5.3) 



with erfc(z) = 1 — erf(z). 

The bound, resonant and scattering state solutions of the potential problem 
characterized by the Hamiltonian Hi can be obtained by solving the Lippmann- 
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Schwinger integral equation. The bound and resonant state wave functions satisfy 
the homogeneous Lippmann-Schwinger equation 

\* l ) = Gf(E)V l \H> l ) (5.4) 

at real negative and complex E energies, respectively. While the wave func- 
tion ) describing a scattering 
Schwinger equation (Section pTg|) 



tion ) describing a scattering process satisfies the inhomogeneous Lippmann- 



I*/ (±) ) = I^ (±) )+G^(£;±i0)^|*r ; ) , (5-5) 



where |^| ± ^) is the solution to the Hamiltonian Hp with scattering asymptotics. In 
Equations (5A), (|5.5|) Gf(z) denotes the radial Coulomb Green's operator defined 



as G?{z) = {z-H?)-\ 

We are going to solve these equations by using a discrete Hilbert space basis 
representation in a unified way by approximating only the potential term Vi. For 
this purpose we write the unit operator in the form 

1 = lim 1 N , (5.6) 

where 

N N 
n=0 n=0 

In this case the {\n), \n)} biorthonormal basis is specified as the Coulomb-Sturmian 



basis ( 4.25 ). The a factors have the properties linin^oo er„ = and lirnjv^oo &n 



1, and render the limiting procedure in ( |5.6D smoother. They were introduced orig- 



inally for improving the convergence properties of truncated trigonometric series 
[IH, but they turned out to be also very efficient in solving integral equations in 
discrete Hilbert space basis representation pS, The choice of 



N _ 1 - exp{-[a(n -N- 1)/(N + l)] 2 } 



° n 1 - exp(-a 2 ) (5 ' 

with a ~ 5 has proved to be appropriate in practical calculations. 
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Let us introduce an approximation of the potential operator 

N 

Vi = 1VJ1 « \ N Va N = V t N = J2 \") (n'\ , (5.9) 

n,n'— 

where the matrix elements 

Z„„'=^«|n>£ , (5.10) 



in general, are to be calculated numerically. This approximation is called separable 



expansion, because the operator V, N , e.g. in coordinate representation, takes the 



form 

N 



(r\V N \r<) = J2 (r\r£) V nn , (n'\r') , (5.11) 
i.e. the dependence on r and r' appears in a separated functional form. 



With this separable potential Eqs. (|5.4| ) and ( p.5| ) are reduced to 



N 



|*,)= Yl Gf(E)\n)V nn , (n'|*,) (5.12) 



and 

N 



\*l ±) ) = \*? ) ) + £ Gp(i? ± i0)|n) Zw (^l*[ ±} ) , (5-13) 

n,n'— 

respectively. To derive equations for the coefficients = (n'\^j ) and = 

(n'\^>i), we have to act with states (n"\ from the left. Then the following ho- 
mogeneous and inhomogeneous algebraic equations are obtained, for bound and 
scattering state problems, respectively: 

[(Gp(i?))- 1 -Z ; &=0 (5.14) 
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and 



where the overlap <& 



[(&(E±iO))- 1 -vM ±) = *t 



(±) 

ill 



n|$| ) can also be calculated analytically 



(5.15) 
I. The 



homogeneous equation ( ]5.14| ) is solvable if and only if 

det[(Gp(-B)) _1 -Zi] = (5.16) 

holds, which is an implicit nonlinear equation for the bound and resonant state 
energies. As far as the scattering states are concerned the solution of (5.15) provides 
the overlap (n\^i). From this quantity any scattering information can be inferred, 
for example the scattering amplitude corresponding to potential Vi is given by jl9| 



•,(+)\ 



(5.17) 



Note that also the Green's matrix of the total Hamiltonian, which is equivalent 
with the complete solution to the physical system, can be constructed as 



G l (z) = [{Gf(z)r 1 -V l 



(5.18) 



Finally, it should also be emphasized, that in this approach only the potential 
operator is approximated, but the asymptotically important H c term remains in- 
tact. The properties of the short-range potential is buried into the numerical values 
of the matrix elements. Thus the method is applicable to all types of potentials, as 
long as we can calculate their matrix elements somehow. Beside usual potentials 
this equally applies to complex, momentum-dependent, non-local, etc. potentials 
relevant to practical problems of atomic, nuclear and particle physics. Further- 
more, the present formalism is equally suited to problems including attractive or 
repulsive long-range Coulomb-like and short-range potentials. The solutions are 
defined on the whole Hilbert space, not only on a finite subspace. The wave func- 
tions are not linear combinations of the basis functions, but rather, as Eqs. ( |5.13| ) 
and ( |5.12| ) indicate, linear combinations of the states Gf(E)\n), which have been 
shown to possess correct Coulomb asymptotics J2£|. 



5.1. MODEL NUCLEAR POTENTIAL CALCULATION 



59 



5.1.1 Bound states 



First we consider only the nuclear part of potential (5T ) and switch off the Coulomb 
interaction by setting Z = 0. According to Ref. |6^|, this potential supports al- 
together four bound states: three with I = and one with I — 2. However, it is 
known that the first two I = and the single I = 2 state are unphysical, since they 
are forbidden because of the Pauli principle. This fact is not taken into account in 
this simple potential model. Although from the physical point of view these Pauli- 
forbidden states have to be dismissed as unphysical, they are legitimate solutions 
to our simple model potential. The proper inclusion of the Pauli principle into 
the model would turn the potential into a non-local one. This problem has been 
considered within the present method in Ref. [E7|. 



In Table |5,1| we present the results of our calculations for the three I = states 
showing the convergence of the method with respect to N, the number of basis 
states used in the expansion. We determined the energies of these states from Eq. 



( 5.16| ), using the CS parameter b = 4 fm 1 . It can be seen that the method is 



very accurate, convergence up to 12 digits can easily been reached. We note that 
according to Ref. |6^| , the energy of the two lowest (i.e. the unphysical) I = states 
is E = -76.903 6145 and E = -29.000 48 MeV in the uncharged case, which is in 
reasonable agreement with our results. 



5.1.2 Resonance states 

Switching on the repulsive Coulomb interaction (Z = 2) the bound states are 
shifted to higher energies. The most spectacular effect is that the third I = 
state, which is located at E — —1.608 740 8214 MeV in the uncharged case, moves 
to positive energies and becomes a resonant state. This is in agreement with the 
observations: the a — a system (i.e. the 8 Be nucleus) does not have a stable ground 
state, rather it decays with a half life of 7 x 10~ 17 sec. 

In our calculations we determined the energies corresponding to this resonance 
and to other ones as well by the same techniques we used before to find bound 
states. In fact, we used the same computer code and the same CS parameter 
(6 = 4 fm -1 ) as we used in the analysis of bound states. The method, again, 
requires locating the poles of the Green's matrix, but not on the real energy axis, 
rather on the complex energy plane. 
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N E 00 (MeV) £iq (MeV) E 20 (MeV) 



8 


-76.903 557 1529 


-29.005 


234 


9134 


-1 


739 478 


2626 


10 


-76.903 609 9717 


-29.000 


352 


3141 


-1 


637 269 


0831 


15 


-76.903 614 3090 


-29.000 


469 


8249 


-1 


608 824 


6403 


18 


-76.903 614 3254 


-29.000 


470 


2338 


-1 


608 742 


5166 


20 


-76.903 614 3263 


-29.000 


470 


2566 


-1 


608 741 


0685 


25 


-76.903 614 3265 


-29.000 


470 


2623 


-1 


608 740 


8256 


28 


-76.903 614 3265 


-29.000 


470 


2625 


-1 


608 740 


8216 


30 


-76.903 614 3265 


-29.000 


470 


2625 


-1 


608 740 


8213 


35 


-76.903 614 3265 


-29.000 


470 


2626 


-1 


608 740 


8214 


40 


-76.903 614 3265 


-29.000 


470 


2626 


-1 


608 740 


8214 



Table 5.1: Convergence of the / = bound state energy eigenvalues E n \ in V a ^ a (r) 
in the uncharged (Z = 0) case. N denotes the number of basis states used in the 
expansion. 



In Table 5^ we demonstrate the convergence of our method with respect to N 
for the lowest I = and I — 2 resonance states. In Fig. 5A we plotted the modulus 
of the determinant of the ( [5.18] ) Green's matrix over the complex energy plane for 
/ = 2. The resonance is located at the pole of this function. Finally, we mention 
that there is a resonance state for I = 4 at -E re s,4 = 11-791 038 — i 1.788 957 MeV. 

Although it is not our aim here to reproduce experimental data with this sim- 
ple potential model, we note that the corresponding experimental values |6^| are 
E teSi i = o =0.09189 MeV, E ies ,i= 2 =3.132 ± 0.030 MeV, E TCS j =4 =11.5 ± 0.3 MeV, 
and r i=0 /2 = (3.4 ± 0.9) x lO" 6 MeV, r /=2 /2=0.750 ± 0.010 MeV, T ;=4 /2 -1.75 
MeV. 



5.1.3 Scattering states 

In order to demonstrate the performance of our approach for scattering states 
we calculated scattering phase shifts 61(E) for V a - a (r) in ( |5.l| ). As described 
previously in this subsection, phase shifts can be extracted from the scattering 
amplitude given in Eq. (5.17). Specifying this formula for the Coulomb-like case 
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N £ rcs , (MeV) £ res ,2 (MeV) 



8 


-0.000 


854 


9596 +i 0.000 


000 


0000 


2.807 


21 


— i 


0.607 11 


10 


0.063 


364 


2503 


-i 0.000 


000 


0681 


2.866 


30 


— i 


0.628 56 


15 


0.091 


785 


0787 


-i 0.000 


002 


8092 


2.889 


68 


— i 


0.620 99 


18 


0.091 


963 


0277 


-i 0.000 


002 


8572 


2.889 


34 


— i 


0.620 53 


20 


0.091 


969 


7296 


-i 0.000 


002 


8588 


2.889 


24 


— i 


0.620 56 


25 


0.091 


971 


8479 


-i 0.000 


002 


8592 


2.889 


23 


— i 


0.620 62 


28 


0.091 


971 


9788 


-i 0.000 


002 


8592 


2.889 


25 


— i 


0.620 62 


30 


0.091 


972 


0064 


-i 0.000 


002 


8592 


2.889 


25 


— i 


0.620 61 


35 


0.091 


972 


0258 


-i 0.000 


002 


8592 


2.889 


24 


— i 


0.620 61 


40 


0.091 


972 


0290 


-i 0.000 


002 


8592 


2.889 


25 


— i 


0.620 61 



Table 5.2: Convergence of the energy eigenvalues E les j for the I = and I = 2 
resonances in the V a - a (r) potential. N denotes the number of basis states used in 
the expansion. 




Figure 5.1: The modulus of the determinant of the G(E) Green's matrix for 
the a — a potential on the complex energy plane for 1 = 2. The pole at E = 
2.889 25 — i 0.620 61 MeV corresponds to a resonance. 
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and for a given partial wave I we have 

ai = i exp(i(2?7; + Si)) sin<5 ; , (5.19) 

where a; is the Coulomb-modified scattering amplitude corresponding to the short- 
range potential, r\i = argr(/ + 1)) is the phase shift of the Coulomb scattering 
with 7 = Z 2 e 2 m/h 2 k and Si is the phase shift due to the short-range potential. 
The convergence of the phase shifts with respect to N is demonstrated in Table 



5.3, where Sq(E) is displayed at three different energy values. As in our calculations 



for the bound and the resonance states, we used b = 4 fm 1 here too. 



In Fig. ^2 we plotted the scattering phase shifts 5i(E) for I — 0, 2 and 4 



up to E = 30 MeV. In all three plots in Fig. 5.2 the location of the corresponding 
resonance is clearly visible as a sharp rise of the phase around the resonance energy 
E Tes . This rise is expected to be more sudden for sharp resonances, and this is, in 
fact, the case here too. The phase changes with an abrupt jump of ir for the sharp 
/ = resonance, while it is slower for the broader I = 2 and I = 4 resonances. We 



also note that the phase shifts plotted in Fig. 15^ are also in accordance with the 
Levinson theorem, which states that <$z(0) = run, where m is the number of bound 
states in the particular angular momentum channel. Indeed, as we have discussed 
earlier, there are two bound states for I = 0, one for I = 2 and none for I = 4. 
As an illustration of the importance of the smoothing factors we show in Fig. 



5.3 the convergence of the phase shift Sq(E) at a specific energy E — 10 MeV with 
and without the smoothing factors in (5.7), and consequently in (5.1C). (Here 



and everywhere else the a parameter of the (5.8) a factors was chosen to be 5.2.) 
Clearly, the convergence is much poorer without the smoothing factors. We note 
that this also applies to other quantities calculated for bound and resonance states. 



5.1.4 Conclusions 

In this section we used the Coulomb-Sturmian representation of the Coulomb 
Green's operator to solve the Lippmann-Schwinger integral equation with a model 
nuclear potential describing the interaction of two a particles. In practice this 
means the approximation of the potential term on a finite subset of this basis. This 
is the only stage where approximations are made (remember the Green's operator 



5.1. MODEL NUCLEAR POTENTIAL CALCULATION 



63 



N £ = 0.1MeV £=lMeV E = 30 MeV 



8 


6.283 


230 


8.817 731 


7.783 217 


10 


9.424 


059 


8.862 581 


4.817 163 


15 


9.424 


018 


8.859 651 


4.835 479 


18 


9.424 


022 


8.859 467 


4.829 861 


20 


9.424 


023 


8.859 441 


4.828 882 


25 


9.424 


024 


8.859 419 


4.828 563 


28 


9.424 


024 


8.859 414 


4.828 555 


30 


9.424 


024 


8.859 412 


4.828 554 


35 


9.424 


024 


8.859 411 


4.828 552 


40 


9.424 


024 


8.859 411 


4.828 552 



Table 5.3: Convergence of the Sq(E) phase shift (in radians) in the V a - a (r) 
potential at three different energies. N denotes the number of basis states used in 
the expansion. 




1=0 
1=2 
1=4 












10 



20 



30 



E(MeV) 



Figure 5.2: Scattering phase shifts Si(E) (in radians) in the a — a potential for 
I = 0, 2 and 4. The resonances in these partial waves appear as sharp rises in the 
corresponding phase shifts. In these calculations a basis with N = 35 was used. 
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Figure 5.3: Convergence of the scattering phase shift Sq(E) (in radians) at E = 10 
MeV calculated with (solid line) and without (dashed line) smoothing factors. 



is calculated analytically), otherwise this method is exact, analytic and provides 
asymptotically correct solutions. Consequently, bound, resonance and scattering 
problems can be treated on an equal footing, while these phenomena are usually 
discussed in rather different ways in conventional quantum mechanical approaches. 
The unified treatment is also reflected by the fact that all the calculations are made 
using the same discrete basis, containing also the same basis and other parameters. 

Finally, we should call the attention upon the fact that this method is very 
accurate. Reasonable accuracy is reached already at relatively small basis, around 
N = 20. The accuracy gained in larger bases is beyond most of the practical re- 
quirements. Test calculations have been performed on a linux PC (Intel PII, 266 
MHz) using double precision arithmetic. The calculation of a typical bound or reso- 
nant state energy requires the evaluation of the potential matrix by Gauss-Laguerre 
quadrature and finding the zeros of the determinant ( 5.16| ) , which incorporates the 
evaluation of the Coulomb Green's matrix and the the calculation of a determi- 
nant by performing an LU decomposition in each steps. The determination of the 



energy value in the first column and last row of Table 5.1, which meant 6 steps in 
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the zero search and handling of 40 x 40 matrices, took 0.06 sec. The corresponding 
resonance energy value in Table required 12 steps in the zero search on the 
complex energy plane and 0.8 sec. The evaluation of the three phase shift values 
in the last row of Table 5.3 took 0.19 sec, 0.11 sec. and 0.08 sec, respectively. So, 
this method is not only accurate but also very fast. 



5.2 An atomic three-body problem 

The bound state problem of the Helium atom is investigated as an example of the 
application of the Faddeev-Merkuriev (FM) integral equations to the atomic three- 
body problems. The FM equations describing three charged particles interacting 
through the long-range Coulomb potential are solved using the Coulomb-Sturmian 
discrete Hilbert space basis representation, which transforms the integral equa- 
tions into a matrix equation. The solution of the matrix equation is possible due 
to the analytic representation of a three-body Green's operator constructed from 
two independent two-body Green's operators by performing a convolution integral 



(Section |2J). 

We recall here that for three-body systems the Faddeev equations H are the fun- 
damental equations. Their homogeneous form is fully equivalent to the Schrodinger 
equation and after one iteration they possess connected kernels, consequently they 
are, in fact, Fredholm integral equations of second kind. Therefore the Fredholm 
alternative applies: at certain energy either the homogeneous or the inhomogeneous 
equations have solutions. Three-body bound states correspond to the solutions of 
the homogeneous Faddeev equations at real energies. 

The Hamiltonian of an interacting three-body system can be written as 

H = H° + v a +v + v~ / , (5.20) 

where H° is the three-body kinetic energy operator and v a denotes the interaction 
in subsystem a (i.e. the pair interaction of particles j3 and 7). Here the usual 
configuration space Jacobi coordinates x a , y a are used, x a is the coordinate of the 
(/3, 7) pair and y a is the coordinate of the particle a relative to the center of mass of 
the pair (f3,j). Therefore v a = v a (x a ) represents the interaction of the pair (/3,7) 
and depends only on the x a relative coordinate of this pair. In Jacobi coordinates 
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the H° three-body kinetic energy operator is given as a sum of two-body kinetic 
energy operators 

H ° = h L + h l = h ° + h l = < + h l • ( 5 - 21 ) 

Suppose the v a , vp,v 7 pair interactions are short-range-type potentials. In this 
case the Faddeev procedure leads to mathematically sound integral equations for 
the three-body system. The three-body wave function |\&) is decomposed as 

|*> = \4> a ) + l#) + |V 7 >. ( 5 - 22 ) 

where the components are defined by 

\i>i) = G°Vi\^) for i = a,/3,7 1 (5.23) 

with the G°(z) = (z — H )^ 1 being the free Green's operator. For bound states 
the set of homogeneous Faddeev equations appear as 

|Va) = G a (E)[v a \1> fi ) + v a \1>J] (5.24) 
1-0/3) = G^v^+vpl^)} 
= G 7 {E)[v 7 \ip a ) + v 7 \^0)} 

where G a (z) = (z — Ha)^ 1 with H a = H° + v a . It is proven that the above set of 
coupled integral equations for short-range potentials have an unique solution. 

Unfortunately in the case of the scattering of three charged particles the situa- 
tion becomes more complicated since the Coulomb potential enters into the game. 
The Faddeev equations originally were derived for short range interactions and if 
we simply plug in a Coulomb-like potential they become singular. The solution has 
been formulated by Faddeev and Merkuriev in a mathematically sound and 
elegant way via integral equations with connected (compact) kernels and configu- 
ration space differential equations with asymptotic boundary conditions. 

In practice, however only this latter version of the theory has been applied, 
since the FM integral equation formulation was too complicated for practical use. 
Moreover, in bound-state problems only the original version of the Faddeev equa- 
tions were applied [Bo, K37L b8| which, in sound mathematical sense, are not well- 
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behaved for Coulomb case and only the bound state nature of the problems helped 
to overcome the difficulties, however slow convergence in partial wave channels was 
reported |6(], |6^] . In order to find a remedy, the equations were solved in total an- 
gular momentum representation, which led to three-dimensional equations So, 
due to these problems and difficulties, the belief spread that the Faddeev equations 
are not well-suited for treating atomic three-body problems and other techniques 
can perform much better, at least for bound states. 

Recently, a novel method was proposed by Papp for treating the three-body 
Coulomb problem via solving the set of Faddeev-Noble and Lippmann-Schwinger 
integral equations in Coulomb-Sturmian discrete Hilbert space basis representa- 
tion. The method was elaborated first for bound-state problems |[l6| with repulsive 
Coulomb plus nuclear potential, then it was extended for analyzing p — d scatter- 
ing at energies below the breakup threshold [Q. In these calculations excellent 
agreements with the results of other well established methods were found and the 
efficiency and the accuracy of the method were demonstrated. This approach has 
also been applied to atomic bound-state problems J6{|. The Coulomb interactions 
were split, a la Noble [[70), into long-range and short-range terms and the Faddeev 
procedure was applied only to the short-range potentials. By studying benchmark 
bound-state problems, contradictory to the conventional approaches, a fast conver- 
gence with respect to angular momentum channels was observed. This approach, 
however, has a limitation which is related to Noble's splitting. The Noble splitting 
does not separate the asymptotic channels, so the equations are applicable only in 
a restricted energy range of low- lying states. 

In this section, instead of the Noble splitting, we make use of the mathematically 
sound Merkuriev splitting of the long-range Coulomb potential [Q, and by doing 
so we solve the Faddeev-Merkuriev integral equations. 



5.2.1 Faddeev-Merkuriev integral equations 

The fundamental equations for the Coulomb three-body problem is discussed here 
following Merkuriev ideas. We recall that the Faddeev procedure is based on the 
observation that the Hamiltonian of the three-body system can be written as the 
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sum of the asymptotically relevant H™' and irrelevant short-range terms 

H = H {1) + v« + v ( ° ] + v« , (5.25) 

where 

HW=H° + vV>+vW+vy. (5.26) 



In equation (5.25) we assumed that the long-range Coulomb potential can be split 



v c = + (5.27) 

by using an appropriate cut-off procedure. 

By applying the Faddeev procedure for the short-range part of the potentials 
we can derive connected kernel integral equations. For this, however it is essential 
that the asymptotically relevant Hamiltonian H"> possesses only continuous spec- 
trum. This property guarantees the asymptotic filtering behavior of the Faddeev 
decomposition |7l), and thus the asymptotic orthogonality of the Faddeev com- 
ponents. A potential with an attractive Coulomb tail has infinitely many bound 
states accumulated at the lower edge of the continuous spectrum, thus the long- 
range part of an attractive Coulomb potential always has infinitely many bound 
states. Fortunately, this is true only in the two-body Hilbert space. In a broader 
space, like the three-body Hilbert space, this statement is not necessarily valid. In 
the three-body Hilbert space, where we have to consider our potentials, we simply 
have an extra kinetic energy term, which can modify the character of the spectrum. 

Merkuriev proposed to split the Coulomb potential (x a ) into a short-range 
and a long-range part by introducing a cut-off function ( a = ( a (x a , Ua) defined on 
the three-body configuration space leading to 

v£\x a ,y a ) = v%(x a )Ca(Xa,y a ), (5.28) 

and 

v®(x a ,y a ) = v%{x a )[l - C a (x a ,y a )]- (5.29) 
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The function £ Q is constructed in such a way, that it separates the asymptotic two- 
body sector fi Q from the rest of the three-body configuration space. On the region 
of il a the splitting function ( a asymptotically tends to 1 and on the complementary 
asymptotic region of the configuration space it tends to 0. Rigorously, fl a is defined 
as a part of the three-body configuration space where the condition 



pal < Uaoi^aO ~l~ \lla 

|/Z/a ) lA, ,with 

*^a0 7 Vol 

o > 0, v > 2, 



(5.30) 



(s) 

is satisfied. So, in fl a the short-range part v a coincides with the original Coulomb- 
like potential v% and in the complementary region vanishes, whereas the opposite 
holds true for v^ 1 ■ A possible functional form for £ is given by 



C{x,y) = 21 1 + exp 



{x/x y 



i + y/ya 



(5.31) 



where the parameters xq and yo are rather arbitrary. Fig. |5.4| shows a typical 
example for the short- and long-range part of an attractive Coulomb potential, 



respectively. Merkuriev proved that H"> of Eq. ( 5.25 ) with long-range potentials 
defined in this way possesses only continuous spectrum if the x a0 and y aQ param- 
eters for all a fragmentations are chosen big enough (see p. 248 of Ref. p8[ ). 

Now, if we follow the Faddeev procedure with the Merkuriev split of ( 5.28|) and 
( 5.29| ) we obtain the Faddeev-Merkuriev integral equations 



|V«) = Gg)(£)[ v «|^) + »«^ 7 )], 



where the G^a channel Green's operator is defined as 



G^{z) = (z-H^)-\ 



(5.32) 



(5.33) 



with 



(0 - ff(') j_„( s ) - 



m> = h 



= H u 



(5.34) 



It is proven that Eqs. ( 5.32) ) below the threshold of the continuous spectrum 
of the Hamiltonian H allow the nontrivial solutions only for discrete set of energy 
corresponding to the binding energies of an atomic three-body system. 
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v (s) 


v (l) 



Figure 5.4: The short- and the long-range part of an attractive Coulomb potential 
in the (x, y) three-body configuration space with parameters xq = 1, yo = 5 and 
v = 2.2 of C defined in Eq. ( p!p ). 



The solution of the Faddeev equations (5.24) and( |5.32| ) necessitates the determi- 
nation of the asymptotically relevant Green's operator in some basis representation. 
This can be done using convolution integral technique of Section 2.4. 



5.2.2 Solution in Coulomb-Sturmian space representation 

The three-body Hilbert space is a direct product of two-body Hilbert spaces, so 
an appropriate discrete basis in the three-body Hilbert space can be constructed 
as the angular momentum coupled direct product of the ( |4.25|) Coulomb-Sturmian 
basis functions as 

\nvl\) a = [\nl) a ®\v\) a ], (n,v = 0,1,2,...), (5.35) 

where I and A denote the angular momenta associated with the coordinates x and 
y, respectively, and the bracket stands for angular momentum coupling. With this 
basis the completeness relation takes the form (with angular momentum summation 
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implicitly included) 



N 



1 = lim V \nvl\) a a (nvl\\= lim 1^. (5.36) 

nM—Q 



It should be noted that in the three-body Hilbert space three equivalent bases 
belonging to fragmentation a, (3 and 7 are possible. 



In Eqs. (5.32) we introduce the approximation 



|V^>=GW(z)1^^1^|V 7 ), (5.37) 

i.e. the short-range potential % in the three-body Hilbert space is taken in the 
separable form 

N 



where v^l = a {nvl\\va\n'v'l'\')p. These matrix elements can be calculated nu- 
merically by making use of the transformation of Jacobi coordinates. The ket and 
bra states are defined for different fragmentation, depending on the environment 
of the potential operators in the equations. 

Now, with this approximation, the solution of the homogeneous Faddeev-Mer- 
kuriev integral equations turns into the solution of a matrix equation for the com- 
ponent vector tp = a {nvl\\ip a ) 

± a = Gg>(z)^J2± y > ( 5 - 39 ) 

where G^) — a {nvl\\G^a \n'v'l'\') a . A unique solution exists if and only if 

det{[G ( ' ) (z)]- 1 ~ 2 ;( s )} = 0. (5.40) 

The solution of the above matrix equation requires the determination of the 
Coulomb-Sturmian representation of the three-body Green's operator Ga of Eq. 
( 5.33| ). Normally this task necessitates the solution of Faddeev-type integral equa- 



2 \rwl\) a 0{n'v'l'\% (5.38) 
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tions or a triad of Lippmann-Schwinger integral equations |[73j. However, the 
Hamiltonian Ha has the peculiar property, that it supports bound state only 
in the subsystem a, and thus there is only one kind of asymptotic channel, the 
a channel. For such a system, if the asymptotics is treated properly, one single 
Lippmann-Schwinger equation is sufficient for an unique solution JT^]. Merkuriev 
proposed an appropriate Lippmann-Schwinger equation which provides a unique 
solution for g£ 

G$(z) = G%{z) + G a J(z)V« s GW(z), (5.41) 

where the three-body potential V£ s is defined such that it decays faster than the 
Coulomb potential in all direction of the three-body configuration space. The 
operators G% s and V£ s are complicated three-body operators, and to construct 
them, Merkuriev used different approximation schemes in different regions of the 
configuration space. 

It is important to realize that in our approach to get the solution only matrix 
elements of between finite number of square integrable CS functions are needed. 
That is the reason why the matrix elements of the channel Green's operator can 
be obtained as 

G®{z) = G a (z) + G a (z)U a G^(z), (5.42) 



where 

G Q = a {r^\G a \r^f\') a , U a = a {nvl\\U a \n'v'V 'A%, (5.43) 

with G a (z) = (z — Ha)^ 1 and U a — Hjp — H a . The channel-distorted long-range 
Hamiltonian H a is defined as 

H a =h c Xa +h«l (5.44) 

where = h° Xa + v% is the two-body Coulomb Hamiltonian. 

Since Ha does not generate rearrangement channels, we should define hy\ 
such that H a also preserve this property. If we are dealing with repulsive Coulomb 
interactions, i.e. e a (ep + e 7 ) > 0, this requirement can be easily fulfilled and hyl 
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is defined by 



n yo, 



h° 



+ e a (ep + e 7 )/y a 



(5.45) 



and U a takes the form 



U a 



e a (ep + e 7 )/y Q 



(5.46) 



On the other hand, in the case of an attractive interaction, like the Helium atom 
where e a (ep + e 7 ) < 0, the attractive Coulomb tail of hy a generate infinitely many 
bound states which has to be pushed away from the spectrum. In practice this is 



achieved by introducing a repulsive Gaussian term into (5.46) 



U° 



(0 



+ vV> - Aexp(-K?/^) - e a (e/3 + e 7 )/y a , 



(5.47) 



where A and k are free parameters. Now the ( |5.43|) matrix elements of U a can easily 
be calculated numerically applying basis functions from the same fragmentations 
on both sides. 

We note here that an alternative and mathematically more sound procedure 
for removing the bound states generated by the attractive Coulomb potential is 
being under development. The basic idea of the new procedure is the application 
of projection operators constructed as convolution integrals of Green's operators. 



The most crucial point in the solution of Eq. ( 5.42 ) is the calculation of the 
matrix elements G a . The operator G a is a resolvent of the sum of two commuting 



Hamiltonians, h^ a and hy a according to ( 5.44 ), which act in different two-body 



Hilbert spaces. As it was discussed in Section 2.4, the three-body Green's operator 
G a equates to a convolution integral of two-body Green's operators, i.e. 



Ga(z) = 



dz ' 9xS z ~ z ') 9yl( z ')i 



(5.48) 



where g^J^z) — (z — ) 1 and gyl(z) = (z — hyl) 1 . The contour C should be 



taken such that it encircles the continuous spectrum of hyl so that g% is analytic on 
the domain encircled by C. For bound state problems this requirement can easily 



be fulfilled (see Fig. (|5.5|)). The matrix elements G a of Eq. (5.43) can be cast into 
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Figure 5.5: Contour C for the convolution integral of Eq. ( 5.48 ). The contour C 
encircles the (continuous) spectrum of gy\ and avoids the singularities of g% . 



the convolution integral of the outer product of the corresponding matrices 
G a (z) = ^r / dz' a (nl\g^Jz- z')\^') aa (^\g2(z')\^') a . 



27ri j c 



(5.49) 



The great advantage of using CS basis is that on this basis the matrix elements of 
the two-body Green's operators are given analytically on the whole complex plane 



(see Section 4.2), thus the convolution integral can be performed in practice. 

As an application, the binding energy of the Helium atom (with an infinitely 
massive core) is determined as the three-body energy to be obtained from equation 



(5.40). The convergence of the binding energy with respect to N, the number of 
basis states in the approximation of Eq. (5.37), is shown in Table [T4|. Excellent 
convergence is achieved up to 6 — 7 significant digits with N ~ 19 applied for n 



and v. In Table p.4| , for comparison, we also quote the result of a very accurate 
variational calculation. 
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N 


I = 


I = 4 


10 


2.897586 


2.903741 


11 


2.897591 


2.903746 


12 


2.897593 


2.903746 


13 


2.897593 


2.903746 


14 


2.897593 


2.903746 


15 


2.897593 


2.903746 


16 


2.897593 


2.903746 


17 


2.897593 


2.903745 


18 


2.897593 


2.903745 


19 


2.897593 


2.903745 



Variational calculation in ref. 74 2.903724376984 



Table 5.4: Convergence of the binding energy of the Helium atom. N denotes the 
maximum value of n and v, i.e. the number of employed basis states. The quoted 
energy values are given in atomic units (h = m e = e 2 = 1). Angular momentum 
states have been taken into account up to / = an d I — 4. The values b — 5, n = 2, 
A = 5 are chosen for the CS basis and the ( 5.47 ) repulsive Gaussian p aram eters 
respectively. The Merkuriev parameters of the splitting function of Eq. ( 5.31 ) are 
xq = 0.7, dq = 15 and v — 2.2 respectively. 



Summary 



In this work we have presented a rather general and easy-to-apply method for 
discrete Hilbert space representation of quantum mechanical Green's operators. 
We have shown that if in some discrete Hilbert space basis representation the 
Hamiltonian takes an infinite symmetric tridiagonal, i.e. Jacobi-matrix form the 
corresponding Green's matrix can be calculated on the whole complex energy plane 
by a continued fraction. The procedure necessitates only the analytic calculation 
of the Hamiltonian matrix elements, which are used to construct the coefficients 
of the continued fraction. This continued fraction representation of the Green's 
operator was shown to be convergent for the bound state energy region. The 
theory of analytic continuation of continued fractions was utilized to extend the 
representation to the whole complex energy plane. The presented method provides 
a simple, easily applicable and analytically correct recipe for calculating discrete 
basis representation of Green's operators. 

The general procedure was applied to determine matrix representation of spe- 
cific Green's operators. The D-dimensional Coulomb Hamiltonian was shown to 
possess a Jacobi-matrix structure on the Coulomb-Sturmian basis. The tridiag- 
onal matrix elements of the Hamiltonian were used to construct the continued 
fraction representation of the Coulomb Green's operator. Numerical tests on the 
convergence of the continued fraction and on the effect of the analytic continuation 
were presented on this example. Continued fraction representation of relativistic 
Green's operators corresponding to the Klein-Gordon and the second order Dirac 
equation could be determined because of the Jacobi-matrix structure of the rela- 
tivistic Hamiltonians on the relativistic Coulomb-Sturmian basis. The relativistic 
energy spectra of the hydrogen-like atoms are calculated as the poles of the rela- 
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tivistic Coulomb Green's operator and demonstrate the high precession numerical 
accuracy of our continued fraction representation. The method has been applied 
for the D-dimensional harmonic oscillator, as well. As a non-trivial example the 
generalized Coulomb potential, which is a member of the exactly solvable Natan- 
zon confluent potential class, is considered. The radial Hamiltonian containing this 
potential exhibits a tridiagonal form with analytically known matrix elements on 
the generalized Coulomb-Sturmian basis, so our procedure is also applicable here. 

Once the representation of the Green's operator in a discrete basis is available 
we can proceed to solve few-body integral equations which provide the real test field 
for our Green's operator in respect of its practical applicability and importance. 

The continued fraction representation of the Coulomb-Sturmian space Coulomb 
Green's operator is used for giving a unified solution of the two-body Lippmann- 
Schwinger equation for the bound, resonant and scattering states. The performance 
of this approach is illustrated by the detailed investigation of a model nuclear 
potential describing the interaction of two a particles. 

As a second application the Faddeev-Merkuriev equations are solved for an 
atomic three-body bound state problem. The solution method requires the evalu- 
ation of a three-body Green's operator, wich is done by performing a convolution 
integral of two-body Coulomb Green's operators. This convolution integral repre- 
sents the real test for our Green's matrices calculated by continued fractions over 
the whole complex plane. 

We can conclude that our general, readily computable and numerically exact 
continued fraction method for determining discrete Hilbert space representation of 
Green's operators is turned out to be valuable in solving few-body problems. 



Osszefoglalas 



(The Hungarian summary of the thesis.) 



Elozmenyek 

A benniinket koriilvevo fizikai vilag mikroszkopikus lefrasat ket kiilonbozo irany- 
bol kiserelhetjiik meg. Soktest vagy terelmeleti modszereket kovetve a sok, il- 
letve vegtelen szabadsagi fokkal rendelkezo fizikai rendszereket mint statisztikus 
sokasagokat kezelhetjiik. Ezzel szemben a nehanytest fizika celja az olyan keves 
szabadsagi fokkal rendelkezo rendszerek vizsgalata, ahol az egymassal kolcsonhato 
objektumok viselkedesenek minel teljesebb fizikai megismereset meg realis celkent 
fogalmazhatjuk meg. A kvantummechanikai nehanytest rendszerek elmeleti tanul- 
manyozasa kozponti szerepet tolt be mind az atom-, a mag- es a reszecskefizika fe- 
jlodeseben, hisz az alapveto termeszeti torvenyek megismeresehez nelkiilozhetetlenek 
a nehanytest modellek. 

A kvantummechanikai nehanytest problema alapveto egyenletei, mint peldaul 
a Lippmann-Schwinger- es a Fagyejev-egyenletek, altalaban integralegyenletek. 
Nehanytest fizikaban az integralegyenletek hasznalata a differencialegyenletek helyett 
tobbek kozt azzal az elonnyel jar, hogy a megoldando egyenletek automatikusan 
tartalmazzak a dinamika szempontjabol alapveto hatarfelteteleket. Eppen ez az oka 
annak, hogy a bonyolult aszimptotikus hatarfeltetelekkel jellemezheto szorasi prob- 
lemak tanulmanyozasa soran az integralegyenletek jelentos elonyt elveznek. Azon- 
ban a legtobb gyakorlati alkalmazas eseteben az integralegyenleteket, azok elonyos 
tulajdonsagai ellenere, mellozik es helyettiik inkabb a megfelelo Schrodinger-egyenlet 
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valamely alkalmas kozellteset hasznaljak. Ennek oka abban keresendo, hogy az in- 
tegralegyenletek a Hamilton-operator helyett annak rezolvenset a Green-operatort 
tartalmazzak. A Green-operator meghatarozasa pedig joval bonyolultabb feladat 
mint a megfelelo Hamilton-operator kozeh'tese. 

Egy adott kvantummechanikai nehanytest rendszer teljes Green-operatoranak 
meghatarozasa ekvivalens a problema teljeskoru megoldasaval, ugyanis a Green- 
operator hordozza a fizikai rendszerrol nyerheto osszes informaciot. Tobbek kozott 
a rendszer energiasajatertekei, kotott, rezonancia- es szorasi allapotai, hullam- 
fiiggvenye, allapotsurusege valamint idobeli fejlodese hatarozhato meg a Green- 
operator valamely reprezentacioja ismereteben. 

A legalapvetobb nehanytest rendszerek, mint pi. a szabad reszecske, a har- 
monikus oszcillator vagy a Coulomb-terben mozgo toltott reszecske, Green-operatoranak 
valamely reprezentacioja az irodalombol mar ismert Jl9| , 

torok szamolasat az alkalmasan valasztott reprezentacio (momentum sajat vek- 
torok, harmonikus oszcillator fuggvenyek, Coulomb-Strum-fiiggvenyek) valamint a 
matematika specialis fiiggvenyeinek bizonyos feltetelek melletti intenziv hasznalata 
tette lehetove. 

A Green-operatorok analitikus ismerete kepezi az alapjat egy, az utobbi idoben 
kifejlesztett kvantummechanikai kozeh'to modszernek is |2^]. A modszer soran a 
Hamilton-operator aszimptotikusan meghatarozo tagjainak (kinetikus energia op- 
erator, hosszu hatotavolsagu kolcsonhatasok) megfelelo Green-operatort egzaktul 
analitikusan szamoljak, es csupan az aszimptotikusan irrelevans rovid hatotavol- 
sagu potencialt kozelitik egy szeparabilis kifejtessel a Hilbert-ter egy csonkolt bazisan 
(a modszer ezert kapta a PSE, Potential Separable Expansion elnevezest). A 
hosszu hatotavolsagu tagok egzakt kezelese biztositja a megoldasok aszimptotiku- 
san korrekt voltat. A modszert sikeresen alkalmaztak a kettest rendszereket lefro 
Lippmann-Schwinger- |2^, |2^, ^|, valamint a haromtest rendszerek alapveto 
egyenleteit jelento Fagyejev-egyenletek megoldasara A modszer alkalmazasa 
soran a Fagyejev-egyenletek eseteben a legnagyobb kihlvast egy haromtest Green- 
operator szamolasa jelenti, melyet kettest Green-operatorok konvolucios integral- 
jakent allitanak elo. 

Mindezek miatt alllthatjuk, hogy valamely rendszer Green-operatoranak anali- 
tikus eloalhtasa a kvantummechanikai nehanytest problemak megoldasa szempon- 
tjabol nagy jelentoseggel bir. 



20, 
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OSSZEFOGLALAS 



Eredmenyeim 

Doktori munkam soran kidolgoztam egy altalanos, egyszeruen megvalosfthato, nu- 
merikusan stabil analitikus modszert a kvantummechanikai Green-operatorok dis- 
zkret Hilbert-terbeli reprezentaciojanak lanctortes eloallltasara. A modszer alapjat 
az a felismeres kepezi, hogy amennyiben a rendszer Hamilton-operatora valamely 
diszkret Hilbert-ter bazison vegtelen szimmetrikus tridiagonalis matrix, azaz Jacobi- 
matrix, akkor a megfelelo Green-matrix egy tetszoleges veges almatrixa a tridi- 
agonalis Hamilton matrixelemekbol es egy lanctortbol szamolhato. A lanctort 
egyiitthatoit szinten a Jacobi-matrixelemekbol nyerhetjiik 

Tovabba a kotottenergia-tartomanyban a Green-operator matrixelemei egy harom- 
tagii rekurzios relacio minimalis megoldasaikent adodnak, ami Pincherle tetele 
alapjan biztositja a Green-operator lanctortes eloallitasanak konvergens voltat ezen 
a tartomanyon. 

A komplex energia-slk szorasi es rezonancia tartomanyaiban a Green-operator 
lanctortes reprezentaciojat a kotott tartomanyon konvergens lanctort analitikus 
folytatasaval adtam meg. Erre azert volt sziikseg, mert ezen energiaertekeknel a 
Green-matrixelemek mar nem minimalis megoldasai a rekurzionak, s igy az eredeti 
lanctort nem, hanem csak annak analitikus folytatasa konvergens. Az analitiku- 
san folytatott lanctorttel a Green-operator teljes komplex sikon vett analitikus 
reprezentaciojat allftottam elo. 

Numerikus szempontbol fontos felismerest jelentett az, hogy az analitikusan 
folytatott lanctort konvergenciajanak sebesseget es numerikus stabilitasat a lanc- 
tortek Bauer-Muir transzformaciojanak tobbszoros alkalmazasa jelentosen javitja. 

Felhasznalva, hogy a megfelelo bazison az alabbi rendszerek Hamilton operatora 
tridiagonalis szerkezetu, a kifejlesztett altalanos modszert alkalmaztam a kovetkezo 
Green-operatorok diszkret bazison torteno lanctortes reprezentaciojanak eloallltasara: 

• D-dimenzios Coulomb-problema Green-operatora a Coulomb-Sturm-bazison, 

• D-dimenzios harmonikus oszcillator Green-operatora a harmonikus oszcillator 
fiiggvenyek bazisan, 

• Altalanositott Coulomb-potencial Green-operatora az altalanosltott Coulom- 
Sturm-bazison, 
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• Relativisztikus Coulomb-Green-operatorok, Klein-Gordon- es a ma- 
sodrendu Dirac-egyenlet, a relativisztikus Coulomb-Sturm-bazison. 

Az analitikusan szamolt lanctortes Green-operatorokat a ket- illetve haromtest 
problemak integralegyenleteinek megoldasara szolgalo kvantummechanikai kozelfto 
modszerekben teszteltem, ezaltal bizonyltva a lanctortes Green-operatorok gyako- 
rlati alkalmazhatosagat: 

• A Coulomb-Green-operator Coulomb-Sturm reprezentaciojanak lanctortes 
eloallitasat hasznaltam fel a Coulomb-Sturm-bazisu PSE modszer alkalmazasa 
soran egy, az alfa reszecskek kolcsonhatasat leiro magfizikai modellpotencialt 
tartalmazo kettest Lippmann-Schwinger-integralegyenlet megoldasara 

• Egy, a Fagyejev-Merkuriev haromtest integralegyenletek megoldasara szol- 
galo kozelito modszerben alkalmaztam a lanctortes Coulomb-Green-operatort 
egy atomfizikai haromtest problema kotott allapoti megoldasa soran. 

Munkam soran a matematikai fizika es a nehanytest kvantummechanika mod- 
szereit alkalmaztam, es eredmenyeim a nehanytest modszerekkel lelrhato atom-, 
mag- es reszecskefizika teriileten alkalmazhatoak. Lehetoseget latok arra is, hogy 
a tridiagonalis matrixok es a lanctortek hasznalatat jelento alapotlet a fizika mas 
teriiletein is sikerrel alkalmazhato legyen. 
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